

Usin9 a linearized theory of thermal-y and r;.ochani- 
cally interacting mixture of linear elastic eolic and 
viscous fluid, we derive a fundamental relation in an 
integral foiiu called a reciprocity relation. Tnia recipro- 
city relation relates the solution of one initial-boundary, 
value problem with a given set of initial and boundary data 
to the solution of a second initial-boundary value problem 
corresponding to a different initial and boundary data for 

a given interacting mixture. Prom this general integral 
* « 

relation we derive reciprocity relations for a heat-con- 
ducting linear elastic solid, and for ^ hear -conducting 
viscous fluid. 

Zn this theory of interacting continua we pose and 
solve an initial-boundary value problem for the mixture of 
linear clastic solid and viscous fluid, we consider tbs 


mixture to occupy a half -space and its motion to be re- 
stricted to one space di^monsion. We proscribe a step 
function temperature on the face of the half -space whore 
face is constrained rigidly against motion. With the 
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aid of the Laplace transform and the contour integration# 
a real integral representation for tho displacement of the 
solid constituent is obtained as one of the principal re* 
suits of this analysis. In addition, early time scries 
expansions of the other field variables are given. 
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INTRODUCTION 


Using a linearized theory of thermally and mechanically 
interacting mixture of, linear elastic solid and viscous 
fluids we derive a fundamental relation in an integral, form 
called a reciprocity relation* This reciprocity relation 
X Jlates the solution of one initial^^boundary value problem 
vith a given set of initial and boundary data to the solu- 
tion of a second initial-boundary value problem corresponding 
io a different initial and boundary data for a given inter- 
dicting mixture. From this general integral relation, we 
derive reciprocity relations for a heat-conducting linear 
elastic solid, and for a heat-conducting viscous fluid. 

In this theory of interacting continua we pose and solve 
an initial-boundary value problem for the mixture of linear 
elastic solid and viscous fluid. We consider the mixture to 
occupy a half-space and its motion to be restricted to one 
space dimension. We prescribe a step function temperature on 
the face of the half-space where the face is constrained 
rigidly against motion. With the aid of the Laplace transform 
and the contour integration, a real integral representation 
for the displacement of the solid constituent is obtained as 
one of the principal results of this analysis. In addition, 
early time series expansions of the other field variables are 
given. 

Chapter I includes a historical survey of early works and 
the various descriptions on mixture theory. 
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Chapter II presents modern mixture theories based on mathe** 
matically sound concepts of continuum mechanics. In Chapter 
III we derive the general integral reciprocity relation for 
a linearized version of an interacting mixture and in 
Chapter XV we pose and solve a basic one«*dimensionaX problem 
using the linearized theory. 



CHAPTER I. HISTORY 

X.l. Darcy's Law 

Tbo th' 20 retical description of the dynamics of 
situationi In which one substance interpenetrates another 
has bc^n a matter of interest to mathematicians, physicists 
and eigineers for many years. The case in Which a fluid 
penrsaces a solid is appropriate to a wide range of problems 
su<^h as soil mechanics, petroleum engineering, water purifi- 
cation, industrial filtration, ceramic engineering, diffu- 
sion problems, absorption of oils by plastics and the re- 
ontry ablation process for spacecraft. A survey of earlier 
worhs on this subject up to 1959 is given by Scheidegger [!]*• 
Ay iarly work on this subject was the study of fluid flow 
through a porous solid with the assumption that the solid is • 
unieformable. Intuitively, “pores" are void spaces which must 
be distributed more or less frequently through the solid if 
t>‘J latter is to be called “porous . “ Extreme small voids in a 
s>lid‘are called “molecular interstices," very large ones are 
c-' lied "caverns." "Pores" are void spaces intermediate be- 
ween caverns and molecular interstices? the limitation of 
:heir size is therefore intuitive and rather indefinite. 

Darcy [2] performed an experiment concerning the flow 
through a homogeneous porous solid. A homogeneous filter 
bed of height h is bounded by horizontal plane areas of 

r — — 

r7urober(s) after name(s) refer to the list of references to 
bo found at the end of this paper. 

This experiment was originally performed by -Darc':j^ in 1856. 
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equal slss A* Theso areas are congruent so that corres«* 
ponding points could be connected by vertical straight linea* 
The filter bed is percolated by an incompressible liquid* 

If open mancmater tubes are attached at the upper and 
lower boundaries of the filter bed, the ligx^id rises to the 
heights and respectively above an arbitrary datum 

level. 3y varying the various quantities involved, one can ' 
deduce the fol3»owing relationship: 

- S<A(h^ - h. ) 

Q ^ (1.1) 

Here, Q is the total volume of fluid percolating in unit 
time, and K is a constant depending on the properties of 
the fluid and of the porous solid. The relationship (1.1) 
is known as Darcy’s law. Darcy’s law can be restated in 
terms of the pressure p and the density r of the liquid. 

At the upper boundary of the bed, whose height is denoted 
by 2 ^ , the pressure is ~ rg (h^ - 22 ^ ' 
lower boundary, whose height is denoted by 2 , , the 
pressure is = rg (h^^ - ) , Here g is the gravitational 

constant. Inserting this statement into (1,1), one obtains 

0 -3 - KA((p2 - P3^)/(rgh) + 1) 

or, upon introduction of a new constant , 

Q =, « k'a(p 2 - ^ rgh)/h * (1.2) 

A constant of the type k' , however, is not very satis- 
factory because one v/ould like to separate the influence of 
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tho porous solid from that of tho liquid. By 1933^ the 
empirical relationship 

k' = k/^ (1.3) 

vas generally accepted whore \x is the viscosity of the 
iluid and H the ’’permeability*’ of the porous solid. 
Physically, permeability measurements are very simple. The 
r?q>erimcnts are performed whereby in a certain system a 
pressure drop and a flow rate are measured. The solution 
of Darcy’s law corresponding to the geometry of the system 
and to the fluid employed is calculated, and a comparison 
between the calculated and the experimentally found results 
immediately yields the only unHnovm quantity k. Darcy's 
law (1.2), when accounting for the separation of the general 
constant into “permeability" and "viscosity," is expressible 
as follov/s: 

g B Q/A = ~ (k/^) {P2 ' Px *** ^9^) A (1.4) 

If the. solid is isotropic and if we consider h as an inf ini* 
tesimal, then the expression (1.4) naturally extends to a 
vector form of Darcy’s law: 

<a = - (k/w) (grad p - rg) (1.5) 

where g is a vector in the direction of gravity. 

Engineering uses of Darcy’s law are limited to flows 
o:diiblting small pressure differentials and to constant 
viscosities and permeabilities. However, for 'liquids at high 
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velocities or for gases, relation (1.1) is no longer valid* 
Further if k and m are variable then this law must be 
modified. 

The validity of Darcy's law has been tested on many 
occasions, and has been shov/n that it is valid for a vidv’^ 
domain of flo*,vs. For liquids, it is valid for arbitrary 
small pressure differentials. It has also been used to 
measure flow rates by determining the pressure drop across 
a fixed porous solid. For liquids at high velocities an<l 
for gases at very low and at very high velocities, Daroy s 
law becomes invalid. 

For given boundary conditions Darcy's law (1.5) is 
itself not sufficient to determine the flow pattern in a 
'^orous solid because it contains three unknowns (q, p, r) » 

1 .V /0 further equations are therefore required for the com^ let© 
specification, of a problem. One is the connection betwe^ \ 
and p of the fluid: 

^ = Jf(p) (1.6) 

e'^id the other a continuity equation, viz.: 

- p If- * div (rq) (1.7) 

where t is the time and P is the porosity defined by the 
fraction of void to the total voliime of the porous solid, h 
great variety of methods for the measurement of the porosity 
are described by Scheidegger [ 1] .The physical conditions of 
flow for which solutions might be sought are (i) steady 

flov/, (ii) gravxty flov/ with a free surface, and (1.- ) 
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unsteady state £Xow« of tliGse# steady state flow solu^ 
tions for incompressible fluids are most easily obtained; 
they are simply represented by solutions of Laplace's 
equation. Except for a few other special cases, Darcy's 
law leads to nonlinear differential equations. 

As an application of Darcy's law we will consider the 
steady state flow of an incompressible fluid* With the 
help of the equations (1.5) and (1.7), one may obtain 

P = div ((rk/n) (grad p - rg) ) (X.8) 

Due to the steady state condition, incompressibility and 
the porous solid being homogeneous, one hasi 

= O (1.9) 

As an example of a steady state solution we give the soiu<» 
tion for two-dimensional radial flow of an Incompressible 
fluid into a well which is completely penetrating the fluid** 
bearing medium. Assuming that the well is a cylinder of 
radiufe , v/ith pressure , and that the pressure at 
distance from the well is p^ , the required solution 

follows easily by considering equations (1.5) and (1.9) as 

2 7T Ic . . 

® log (Rj^/Rq) 

'where Q is the total discharge per unit time. 

A major limitation in this theory is due to the 
assumption that the solid is rigid, in most applications 



8 


this is simply not true. incorporate the effects that 
a deformable solid imposes upon the flow, it is necessary 
to develop some connection between the stresses and the 
corresponding strains of both the fluid and the solid. We 
v/ill consider the case of soil consolidation. 

4 

1.2. Biot’s Work 

A soil under load does not assume an instantaneous 
deflection under that load, but settles gradually at a 
variable rate according to the load variation as in clays 
and sands saturated with water. A simple mechanism to 
explain this phenomenon was proposed by Tersaghi [3] by assum* 
ing that the grains constituting the soil are bound together 
by molecular forces and constitute a porous material with 
elastic propex'ties whl3,s the voids of the elastic skeleton . 
are filled with water. A load applied to this system will 
pr*^^uce a gradual settlement, depending on the rate at which 
tho water is being squeesed out of the voids. Terzaghl 
pppt'ied these concepts to the analysis of the settlement of a 
colu n of soil under a constant load and prevented from 
late al expansion. Tiie remarkable success of this theory in 
pred: 7ting the settlement for many types of soils has led to 
the i c tens ion to the three-dimensional case and the establish*^ 
ment '•f equations valid for an arbitrary load variable with 
time. will review extensive work done by Maurice A* Biot 
in th:*'^ field. 



Biot [4] assumed the following basic properties of the 
soils <1) isotropy of the material# <2) reversibility of 
stress-strain relations under final equilibrium conditions# ' 
(3) linearity of stress«-strain relations# (4) small strains# 
(5) the water contained in the pores is incompressible# (6) 
the water may contain air bubbles* (7) the water flows 
through the porous skeleton according to Darcy's law* We 
refer the points in this continuous medium to a rectangular 
c artesian system# #■ i » 1#2,3. Consider a small cubic 
element of the consolidating soil# its sides being parallel 
the coordinate axes* *2^is element is taken to be larga 
or TUgh compared to the size of the pores so that it may be 
treated as homogeneous#' and at the same time small enough 
•■jropared to the scale of the macroscopic phenomena in which 
we are interested# so that it may be considered as Infini- 
tesimal in the mathematical treatment. Physically the 
stresses of the soil are composed of two parts; one which is 
caused by the hydrostatic pressure of the water filling the 
pores# the other caused by the average stress in the skeleton 
They must satisfy the well-known equilibrium conditions of o 
stress field, Let denote the stress components and lot 

X. denote axes of the cartesian system. 

By we shall mean the jth stress component of the 

skeleton acting on the face x^ constant* Then according to 
the equilibrium for the infinitesimal element of volume we 
have 
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and 


a. . . 

3*30 


O 


(1-iO) 



a.ii) 


Denoting by the component of the displacement in the 

and assuming the strain to be small, the 
values of the strain components are 



1 

2 




( 1 . 12 ) 


In order to describe completely the macroscopic 
condition of the soil, an additional variable giving the 
amount of water in the pores is considered. The increment 
of water voXxamo per unit volume of soil is called the varia* 
tion in water content and is denoted by e , and the incre- 
ment of ^ater pressure is denoted by a . Let us consider a. 
cubic rlement of soil. The v^^ater pressure in the pores may 
be considered as uniform throughout, provided either the 


^ize of the element is small enough or, if this is not the 
case, provided the changes occur at sufficiently slow rate 
to ren-"-r the pressure differences negligible. Since it iv 
that the changes in the soil occur by reversible 


processes, the macroscopic condition of the soil must be a 
definite function of the stresses and the water pressure, 

, l.e.j the seven variables , 0 must be definite funp- 

tions of the variables and a , Furthermore if the 


All s'^h^cripts run over values 1,2,3, and, when repeated 
indicr-o 0 svm on the index over 1,2,3. The notation f. 


vari^hi respect to the jth independent 


variable, i.e., df^dx. , 

3 * 3 



strains and the variations in water content are assumed to 
he small quantities, the relation between txio sets of 
variables nay be taken as linear. Consider the case Where 
a “ o. Tho sl3^ components of strain are then functions only 
of the si:: stress components J. . * Assuming the soil to 
have i":otropic properties, these relations reduce to the 
amII kno^rn escpressions of Hookes' law for an isotropic 
olastio body in the theory of elasticity 


1 , V ft % 

2Q ^^ij ” (l-fv) ij^ 


( 1 . 13 ) 


There the constants G, v may be interpreted, respectivoly^ 
as the shear modulus and Poisson's ratio for the solid 
skeleton. 


The effect of the water pressure a is now introduced. 
By reason of the assumed isotropy of the soil, this effect is 
limited to a dependence upon the three strain components 


^11' ^ 22 * ^33 such dependence is uniform in each direction 

Hence taking into account the influence of a , the relationr 

(1.13) become 




(l+v) 



• ( 1 .: 1 ) 


where H is an additional physical constant which plays the 
rolo of a bulk modulus. These relations express the six strain 
.components of the soil as a function of the stresses in the 
soil and the pressure of the i^ater in the pores . 


•frit 

The Kronecker delta, denoted by 6. . 

..... 3.3 

xt 2. « 3 and O if i j , 


is defined as 1 


X2 


derive the dependence of the increment of water 
content 6 on these same variables Blot considers the 
^«neral relation 


(1*15) 




®3®33 ■*■ °4''l2 “5*^23 * ®6®13 * ' 


and argues that because of the isotropy of the material a 
.•hange in sign of ^^.2' ^23* ^13 affect the water 


content. Therefore 


a. » O and the effect of 
6 


shear stress components on 6 vanishes* Furthermore r 
;'.l three directions must have equivalent 

properties so that =* . Relation (1*15) may be 

written In the form 

-• X 


vrhere are R are two new physical constants. 

To this point in the derivation Biot has used assump^ 
tions (1) ^ (3), (4). He now uses (2) to show that the five 
constants can be reduced to four. This assumption# i.e., 
the existence of a potential energy# means that the work 
?one to bring the soil from the initial state to its final 
3tate of strain and water content is independent of the way 
by which the final state is reached and is a definite 
function of ^ only. The potential energy of the. 

*r 

per unit volume is 


U 


2 


+ 00 ) 


a. 17 ) 
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As a result of seme. elementary manipulations, Biot shows 
that H » , and we may write the equation (1«16) as 

Halations (1*14) and (1.16) are the fundamental 
relations describing completely in first approximation the 
properties of the soil, for strain and water content, under 
equilibrium conditions. They contain four distinct physical 
constants Of v, H and R. Solving equation (1.1^) with 
respect to the stresses, then substituting into the equilib« 
rium conditions (l.lo), one obtains 


.with 


G7 + 


G Se 

1-2 V dx- 
^ 2 . 



0 


<1.19) 


a 


G 

3(l-2v) H 


( 1 . 20 ) 


Th^re are three equations with four unknown u . , a . In 
order to have a complete system, one more equation is needed, 
This equation is derived from Darcy’s law governing the flow 
of water in a porous medium. An elementary cube of soil is 
cons Idered and the volume of water flowing per second per 
unit area through the face of the cube perpendicular to the 
x^-axis is denoted by . According to Darcy 'e law these 
three components of the rate of flow are related to the water 
pressure by the relations 


V 


1 


k ha 


a >21) 



whoKo t’:'" physisaX constant k is the coofSicienfc of 
pemea^ility of the soil, and ji if? the viscosity of the 

Ginco the water is asoumed to be Incosnpreosible# one 


obtains 


M 

dt 




(X 




Fron equations (1*17)# (1.20) and (1.21) on© obtains that 


where 


]S 




a a 



1 ^ 
Q ^t 



(1.23) 


four differential equations (1<,18) and (1.22) are the 
basic equations satisfied by the four unknown u^, o « 

In a paper by Biot and Willis [5], methods of m©asure<- 
ment for the four distinct physical constants v, H and , 
JX are described and the physical interpretation of the cOTi- 
stants in various alterf::atQ forms is also discussed. 

'In a later i^ork by Biot [6] tiia stress-strain relations 
which are valid for the cas© of an elastic porous medium with 
nonunifom porosity, i.e., for which the porosity varies frtm 
point to point are derived and these relations lead to the 
siv equations for the si?s components of the unknown disolac©-*- . 
ra vector fields ta for solid component, 5 for fluid 
rc •'■^onpnt , Tae stress-strain relations are 



an 6ij (Xa 


OK5 ) 


cs 


(l*2^a) 
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« aMe HI (1«24|) 

whore the coetticients M« X> h, t and of equatlono 
(l«24a) and (l«24b) are e^lvalent to the constants *0# 
2Gv/(l«v) , G« 6 and a of the equations of (X«19) and 
(1.23)* the el)c equations for the six components of the 
unknown vector fields u and £ are 

9 

2 (n <Xe - aM?) '• 0 

<k/u) grad (aMe - M5) * 0 

Wo will consider a uni form porosity case, i«e«f let us con*« 
sider a particular case where the coefficients t) » X# a» 

M and k/u are constants. Xn this case equations (1*2S) 
become 

2 

ri V a + (tI ♦ X) grad e - a M grad ? • 0 {l«9da) 

^ » (k/M)Ma grad e - (kM/u) grad ? . (l«2|b} 

these equations can be written in the form of equations 
(1.15) and (1.23) by the application ef the divergence operaWer 
to the equation (1.26b). With the aid of the general Papkovitlio* 
Boutsinesq solution for Lane’s equations of the theory of 
elasticity, the general solutions for the equations (1.26) are* * 
obtained as 

r grad (1^ + r. \ «“•* ♦ U-”®) 

C • trad 4 - grad «iv fj^dt 


{l.iTb) 



16 


Whe::w and are aolutlona of Laplace's equation* 

r l8 the pv^^Mon vector and satiaflee the diff« lion 
equation 

1 Oc/n) . 

Utilizing equations (1.24a) and <l«24b)# and cone, der* 
ing the dynamical case* Biot [7] established equations for 
acoustic propagation in the elastic isotropic porous rr\id 
containing a viscous fluid by adding suitable inertia t^rnis * 
in the previous theory* and discussed the propagation o^ 
three kinds of body wavea. For simplicity of notation % i 
will use a new set of coefficients which are related wi^h the 
coefficients of equations (1.24) i 

H a T1. A • X + B > P(a-P)M. C a (1.28) 

Where P denotes porosity, with the vector notation 


u w (u„* 
B a 


X' \ 


* Ujg) 




NV^u + grad( (A’i'N)e + Be) * 


ht 


xr 


grad (Be + Cc) a — j ^^12^* * ^22®^ 

&t 


Mb 29a) 
M.29b) 


where div u « e, div U « « , and laasa 

coefficients which account for the fact that the relative 
fluid flow through the pores is not uniform. Applying the 
divergence operation to equations (1.29)* one obtains 



X7 


7^<(A+2N)o + B«) « (1.30a) 

V^(Ba ■!• Ce) - (Pj^ 2 « + P 22 «> (1.30b) 

These two equations govern the propagation of dilatational 
waves which involve coupled motion in the fluid and the 
solid* Similarly, applying the curl eperation to equationa 
(X.24) one obtains 

- NV» (1.31a) 

j2 

^ (Pl2« + P 22 O) - O (1.31b) 

Where 


curl u w u> , curl u « 0 . 


These equations govern the propagation of pure rotational 

waves. But there is only one type of rotational wave becaus# 
equations (le31) reduce to 


N7^« a Pll<l 


-fi2 . a^u> 

‘’ii‘*22 at^ 


a 



U) 


(1.32a) 

(1.32b) 


An additional result found in [7] is that there is possibly 
^ a wave such that no relative motion occurs between the fluid 
and solid when a certain relation is satisfied between the 
elastic and dynamic constants. 
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‘ Another type of ral^sture theory was considered, by 
Truesdell and Toupin [8]* In thie theory the concept of 
superimposed continua is introduced? i.e., it is assumed 
that the neighborhood of each point of a material is 
occupied by all membero of the mixture. Xfe define the den- 
sity of the mixture to be the eum of the individual denaitleo 
of each constituent. The velocity of the mixture is defined 
by the requirement that th© mass flow of the mixture is the 
sum of the individual mass flows. Then the position of each 
particle of the mixture is defined by an integration of the 
velocity of the mixture? but such particles, in general, bear 
no sample relation to the particles of the constituents. 

The main results of the work of Truesdell and Toupin are 
the following • 

(a) The mass of the mixture satisfies an equation of conti- 
nuity if the mass cupply of the mixture is ^ero. This equa- 
tion of continuity is precisely that found for ordinary 
continuum mechanics. 

(b) Let the total stress of the mixture be defined as the 
sum of partial stresses plus the stresses arising from dif- 
fusion. Then a necessary and sufficient condition that 
Cauchy's first law holds for th© mixture is t^at momentum 
suppliocl by unbAlanccd inertial forces of the several con- 
stituents plus momentum supplied through the creation of 
constituent diffusing masses oha?!.l add up to zero. 

(c) We define tho Internal energy cf the mixture is the cum 
of the internal energies of the conotituento pluo tho kinetic 
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energies of diffusion. This definition leads to the fact 
* 

that the energy supplied by an excess internal energy rate# 
plus tho energy supplied by the work of the excess inertial 
forces against diffusion, plus the energy due to mass supply# 
i.ust add up to zero for the mixture. 

TruQsdell and Toupin*» work on the mixture theory was 
incomplete in spite of the above results. However, their 
work inspired a nunOber of researchers who have since made 
these theories more complete. For example, Adkins [9}# [10]# 
[11], Green and Adkins [12] among others have given discus- 
sions concerning nonlinear constitutive equations. Kelly [13] 
has extended this work to include electromagnetic effects 
while others have accounted for chemically reacting mixtures. 

Recently, conceptually more simplified theories have been 
developed. The basic equations of mass and momenttom balance 

a 

in these theories are equ 2 .valent to those proposed by Tansesdell 
and Toupin. These theories have led to the formulation of 
linearized equations governing thermome chan leal disturbances. 
Since the present work is within the framework of these 
theories, these theories will be reviewed. 



CHAP^En SX. •:?HEORY OF SNTERACWNG CONTIKUA 
2*1* ITonlin'aar C?heo'*y 

For simplicity, attention ia confined to two constit-* 
uent continva in the theory [14]. T?a consider a mixture 
of two continue and s^ which are in relative motion 

to each other. We will agree to call . a solid and »2 
a fluid. We assume that each point within the mixture is 
occupied simultaneously by a^ and s^# and refer the 
motion of the continua to a fixed system of rectangular 
cartesian axes. The position of a typical particle of 
at time t is denoted by x. (t)# where 

3 . ♦ 

V 

t) (^«<T<t), (2*1) 

X. is a reference position of the particle, and lower and 

A 

upper case Latin indices take the values 1,2# 3. We use 
the notation 

^i “ (2.2> 

and can express (2.1) in the alternative form 

x^(t) » x^(x^,X2»x^, t# t) {2*3) 

i5x. (t)j (t) 

where 1 .^ ^ >0. (2.4) 

I 

Similarly# for a typical particle of s^# we have 


20 
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t), ^ ("“<▼;£*) ( 2 .s) 

or 


y^(‘r) = yi<yi*y2'^3'*'"^ 

together with 


Sy^ C t) f 

|8yi(T) 


\iyj 


( 2 . 6 ) 


(2,7) 


We assume that the particles under consideration occupy 
the same position at time t so that 


y . X. 


( 2 . 8 ) 


Velocity vectors at the point zz. « y . in s. and s^ 

& A «L 2 

at time t are 


( 1 ) 


D X. 

«i “ dF^ ' 


( 2 ) 

0 Vi 

''i “ dT“ 


(2.9) 


( 1 ) 

where /Dt denotes differention with respect to t 

holding fixed in continuum s- and D^^^/Dt denotes 

# J . 

a similar operator for s^# holding fixed. These 

operators may also be written as 

\ 


Dt 


cs 


■ ■^ 11 ju 


u 


m 


m 



&t 


+ V 


h 


( 2 . 10 ) 


Acceleration vectors at time t 

A 

g^. where 



are denoted by 


and 
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(2. XX) , 


The densities of and s^ at tine t are» respec- 
tively# and p^# and the rate of deformation tensors 

at tine t ' are defined to he 


2<*ij " «i. j ^ “j.i ' 


j ♦ ''j.i 


where a comma denotes partial differentiation with respect 
to or yj^. We also define a mean velocity w^ be 

the equation 


W. " Pj_U. P2V. , . P - Pj_ f P2 


(2.13) 


and put 


Dt * It 

m 


(2.14) * 


It then folios that 


^1 Dt “ P ru^ 


( 1 ) 


2 Dt 


( 2 ) 


D 
Dt • 


(2.15) 


Let SB be an arbitrary fixed closed surface enclosing a 
volume B and let n^^ be the outward unit normal to 
Let U be the internal energy per unit mass of the mix- 
ture. The externally applied body forces per unit masses 
of Sj^ and Sj aJf® denoted# respectively# by the vectors 
F . and G. • And these vectors are defined through their 

X X 
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rate of work contributions arbitrary 

velocity fields arid . The surface force vector 

t^ per unit area of 3B is such that the scalar 't^u^ , 
for arbitrary , is a rate of work per unit area of dB • 
And a similar definition can be made for tlie vector 
associated with the velocity vector v. • The scalar r 

3L 

is the heat supply function per unit mass of the two cdri- 
tlnua due to radiation from the external world and heat 
sources. The flux of heat across dB is denoted by a scalar 
h per unit area and unit time. 

Uheorem (Green and Nachdi) 

Let us postulate an energy balance at time t In the 

form 

K 

“ J + P2®i'^i>'^ 

V 

+ *<■ - J’ h dA (2.16) 

A A 

Then itSfollows that . 

it P'^k.H “0 
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which states that maso olemontG of the mixtures aro con- 
served. vectors t^, are defined with referenco to 

an arbitrary s'lrfaco A. \fhen the si'»rface at a point 
is perpendicular to tho Xj^-axis, wo denote the correspond- 
ing valu^^ by and refer to those as stresses* 

Hien it also follows from (2.16) that 

<\i * ’^kt> 'k * Pl^i ■*• ^2®i “ 


P2’'iV 

K 


Which is the equation of motion. 

How let the mixture be composed of an elastic solid and 
‘ a non-Newtonian viscous fluid, •she Kinematic quantities 
entering into the theory for the solid are the velocity n$ a 
strain tensor e, a rate of deformation tensor d and a 
vorticity tensor F; and, for ths fluid, the velocity v# 
a rate of deformation tensor f and a vorticity tensor 
A . 'In addition to the body forces previously defined we 
have for the solid constituent tlie following mechanical 
quantities: a partial stress tensor O for the solid 

and a similar tensor ‘Jr for the fluid. Due to the inter- 
action of the two constituents, the theory gives rise to 
a diffusive resistance vector ^3 • ^She thermodynamic 

# 

quantities, referring to the mixture as a whole# are the 
temperature T# the specific entropy S, the specific 
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Helmholtz free energy A, the heat ^;lux vector <2 and the 
heat supply function r. i?or such a mixture the constitu- 
tive equationo given [15] in carcnical form as 




A « A(e, T) ■ 

(2.19a) 

S 8(e, %') 

(2.19b) 

q « q(®, grad 1 .) 

(2.19c) 

® » u> (grad ", grad 0 , 0 ^, T) 


•f- 0 ( 0 , f, r-A, tt-v, p . T) 

(2.19d) 

cr « 0 ( 0 , ' . f, r-A , u-v > T) 

(2.19«) 

7T « 7r(o, a, f, T-A , u-v, p,, T) 

( 2 . 19 £) 


Farther restrictions upon equations (2.19) arise from a 
•general principle of invariance under superposed rigid 
body motions and a material invariance associated with the 
assumed isotropy of the solid constituent [16,17]. It has 
been shown that in order to arrive at a determinate linear- 
ized theory it is svifficient to adjoin to the linearized 
forms of the field equations (2.18) a system of linearized 
forma Of the constitutive equations (2.19). We nov/ examine 
the linearized theory more closely. 

2.2. Linearised ^eory 

^ Assume that the mixture undergoes a disturbance in 
which:* (a) the material points of the solid constituent 
*are displaced by only small amounts from their positions 
in an equilibrium state of the mixture in which the densi- 
ties of the solid and fluid constituents and the tempera- 
ture have the uniform values respective! 



and (b) the epeetl o^; the fluid constituent io small. We 
refer tho points uith respect to a fixed system of rec- 
tangular cartesian ccxsrdinatoo . 

We consider a thermodynamic process in which the 
motions of the solid and fluid constituents of the mixture 
and the temperature fie. Id ^ each admit pov/er series 
represertatioofi in terms of a positive real nuittoer We 

choose c to be a measure of the extent to which the mix- 
departs from some reference state* As our reference 
conf'p'iration t/e take the eoui librium state of the mixture 
in w"'ich 

^ rs K » Vfi T*=7 (2«20) 

t’* 

•jid 

S* « G = O, r c O (2 •21) 

If a(x#t) X3 a field quantity^ wa denote its €-expansion 
by s^(X,t); thus 

« 2 e” (2.22) 

We assume the following G*~oxpansions 

\ 

X =■ , y = + T = T+ e (Jc.t). (2.23) 

From equations (2.9) e"expansions for the velocity vectors 


u are 
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n « (2^24) 

We suppose that each of the expansions (2«23) and (2«24) 
is absolutely convergent with an interval of convergence 
O < ^ The linearization of equations \2*9) to 

(2*15) is obtained by replacing each of the variables by 
the ^-expansions given in (2.22) to (2*24) • Without 
further details it follows that the linearized dioplacetneht* 
strain relations for the solid are 

5w . ^ 

• <2.2S) 


Prom (2.17) 


P 


1 


the individual continuity equations become 

^2^pp - <2.26) 


The vorticity components of the two constituents are 

given in terms of the velocity components by 


r, ^ Uf , A;.; ® , 

ij li#3J 13 li#j] 


(2.27) 


The application of the principle of invariance under 

superposed rigid body motions of the mixture to equation 

(2.18) leads to the follcifing equations of motion; 

\ 

V.p ’ “i ’"pi.p <^2®i “ ^2^ <2*28> 

On entering the ^r-cxpansions for the various field 


*We use the notation A,, » J(A. ^+A. .), Af. 

Up 3) ** i»3 3#i [ipjj 

»<*i. j • 


28 


quantities Into oquationo (2.16), (2.JL7), (2.26) and 
(2.23) an'1 equating tho coefificients o2 6 , we obtain the 
linearized field equations for the mixture as follows. 


Equations of motions 


a . 

^l^i ** ^1 


in 

ft 

B®, S 2- 0 

(2*29) 

IT 4 + 

pi#p 

«»i P2®i “ ^2 

ht 

in 

B®, t >0 

(2.30) 


Energy equations 


. TT ST . , — 




in B , t > 0* 


(2.31) 


In (2.31), p *o -5- p^ is the total initial density of 
the mixture and w. . ir. * are the diffusive resistance 

* *3 13 

and the partial stress tensors in the equilibrium state. 

The liftearized constitutive equations obtained from (2.19) 
are 

pA « p A + ai«pp + 03<P2 “aCr-T) (2.32) 

pS B- ( 03 + OgSpp ®10<P2"P2^ + O^(T-T)) (2.33) 

“ -J«T.£ - K* («^ " Vj^) (2.34) 

**^*^****~TTerrTi — ITT — TM— rwiM-mrTrw^~nrr— — FT)>ir~>n 11 I iSM ins^iiiinie irriir—^ i iiiaiwi m— 

*Se» notations used in Theorem on Page 31. 
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i' 


2 




p 


■» 9 •*» <*♦ tj» 


+ a"- (^ *• A )* 


(2«35) 




•(ij) “ <«i - <f -“4^ V^‘= ♦®e> <i>2 “ Pa)-»'09<*‘'»> 

•*• + h. V ®i j ♦ (°i •*■ <»5> «i j 


+ aujfij 


( - p ,o, + P, Oj - ag) 6 - 

p p 


(2.36) 


P+P. 


(i j) " ' ■ ‘'2“2 ^2 ^*2 ~ "e' '^pp - '^F~ ®2 * <’2®6’ <<*2"'*2> 


- P2«10 <’■"’’ ^ ^4**pp + j * 2^4di j+2itfij (2.37) 


®fij] “■’^[ij] ” ‘’^ijp%“V "%j* 


(2.38) 


We note that there are total of 24 constants 
\l# X^'OCa# M# y^# u^? a#a' * pD,D* "pkpK® which have to be 

determined by an experiment for the mixture. The entrOpy*- 
production inequality [14] imposes restrictions upon the 
constitutive equations i^ihich# in the linearized theory# 
require the material constants to satisfy the following 
inequalities : 


*We use the notation *“ '^[ijj 



3C 


■{3K3+2m3+3X4 + 2u^) 2 < 4(3X3L + 2|ij)(3X-!-2p) (2.39) 

a>0, D** >0,(a" -D)^ <4aD",k >0,K’^ <4T«k 

This Xinearized theory is well posed in the sense that the 
number of the field and constitutive equations equa.ls the 
number of field quantities to be determined. We would ex- 
pect that the boundary conditions for the initial-boundary 
value problems for the mixture are similar to the classical 
boundary conditions for the elasticity problem in which 
the stresses, strain, and displacements are sought* We 
\ recall that the classical boundary conditions are: (aj the 

forces may be given on the surface of the body, (b) the 
displacements may be given on the surface of the body, 

(c) the forces may be given on some portions of the body 
surface, while the displacements are given on the other 
portions. Indeed, the proper form of the initial and 
boundary conditions v/hich should be adjoined to the field 
and constitutive equations of the linearized theory of 
interacting continue so that the sufficiently smooth solu- 
* tlons of the field and constitutive equations are deter- 
mined uniquely are quite similar to the classical boundary 
conditions of the elasticity except that we have the 
temperature terms and we have to specify the boundary con- 
ditions to each component of the mixture. These conditions 
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are specified in the following theorem. 


Iheorcin (Afkin. ChadwlcH, Steel) [ 18 ] 

Let a be a bounded regular region of three-dimensional 
Euclidean space occupied by a mixture of an elastic solid 
and a viscous fluid undergoing a disturbance of small 
amplitude during the time interval t ^ O. We denote by 
ZB the boundary and by the Interior of B. We use 

notation ZBy ZB^ and ZB^, ZB^ for arbitrary subsets 
of ZB and their complements with respect to and n 

refers to the unit outward normal vector field on • 
Suppose that the constants « 

k and k' satisfy the conditions (2.39) and that 

and satisfy the inequalities 

\ 


o ^^2 1 2 

Oj^+Og ^ O' f «2'^“6 ^ °'®7 ^ "3> ®4 3®5 

^ ^ a + Og) 2 ^ (a - 5)+04 +|o5) (=05 * “e’ 

P ^P P P 


io. 

( 2 . 40 ) 


Then there exists at most one set of functions pj of 

class and w., T of class Which satisfy equations 

(2^.-25), (2.26), (2.27), (2.29), (2.30), (2.31), (2.34) 
to (2.38) and the subsidiary conditions 

Wi « u^ o v^ = ^2 “ ^2 ^2* ^ 


on B at t s o # 


(2s41) 



32 


Ui V. - E. , (Cp, + ,Tp.)np - on , 




(2.42) 


T a T + O on SB,, < 2 i^n « F on SF, for t ^ O , (2*43) 

<6 P P ^ 

A A A A A 

where w , , u . , v. , p^* ", R. r ^4 # 0^ # V . , © , P and F, , 

4-l.S.X X 

G^ , r arc prescribed functions on the appropriate 
domains and p^, p^* are given, ^strictly positive, 
constants. 

It is i/eli Known that for the dynamical motions of 
linear isotropic elastic solid, the displacement vector w 
may be represen'fed as a siun of two components representing 
motions of dilatational and rotational types, i.e., 

w = grad cp + cvirl ^ , 


where cp and ^ satisfy the wave equations in xdiich 
appear the speeds of propagation of dilatational and rota- 
tional body waves respectively. Ihis representation is 
known to be complete in the sense that every sufficiently 

smooth solution v/ of the equation of motion of linear 

\ 

isotropic elastic solid is expressible in the stated form 
where the scalar and vector functions ^ satisfy the 
above mentioned wave equations and in addition, div ^ ^ o* 
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For the motions of the interacting continua of an 
elastic solid and a viscous fluid, Atkin [18] has estab«* 
lished a decomposition of its motions into components 
representing motions of dilational and rotational typesf« 
Each part of tho decomposition is somei^hat simpler in 
form than the original system of the differential equations 
(2.26), (2.29) to (2.31) and (2.35) to (2.38), but con- 
siderably more complicated than the wave equation. The 
main merit of this new formulation for the motions of the 
interacting continua is that it allows the investigation 
of the propagation of small amplitude plane waves in a 
non-heat conducting mixture of an isotropic solid and an 
inviscid fluid V?e define new material constants by the 
following combinations of the material constants. 


H P 

Kj ® ZZ ^ ^1 * ^1 ^ *44) 

P P 


Introducing the vector differertial operator 

*=* § 9 ^^ad div - n curl^ (2.45) 

and supposing that body forces' and heat sources are absent, 
the governing equations take the form 
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•mm « 

p2 + p2 V a O 

Wi + JD“]w + J^U3+2 , U3-iD“]v 

* M 

- (a + p curl) (w - v) + Hsf + j 
-grad(Kjp 2 /^ + Bj^e/T) f» p ^ w 

3^[X4+2^4, P 4 - ^D«]w + R[X 2 + 2 M 2 ^ ^i 2 
+ <a +^ curl) (w - v) +l.[K^, 0 ]w 
-grad(K2 p2/P2 + B2®/^^ ® ^2^ 

6 + div(B^w + B 2 V - K* (v - v) ) « k 0 

Equations (2.46) contain twenty material constants of which 
nine# the \*s# u*s and D*\ may be regarded as viscosity 
coefficients# the three K*s as bulk moduli, G, as a 

JL 

shear modulus# and B 2 as products of bulk moduli# 

Of the remaining five constants, K* is associated with 
the transfer of heat in the mixture due to the relative mo- 
tion of its constituents, k is the thermal conductivity, 

I 

c^ is the specific heat at constant deformation# a# 0 
, arise from the interaction of the two constituents through 
the diffusive resistance and antisymmetric parts of the 
partial stress tensors. 

— — 

A dot above a variable means partial derivative of the 
variable with respect to time. 


(2.46e) 

(2.46d) 


(2.46a) 


(2.46b) 
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Suppose nc^ that there are ocalar functions p^(x»t) and 
q>2(Xft) which satisfy the equations 

(Xj_ + 2iij)7^9j + (Kj + j + (Xj ♦2 u^V*<j)2 -o(^« 


-KjCPj/p’j - 1) » O (2.47«) 

< X 4 + 2 U 4 > + Kj 7^ <Pj ♦ ( Xj+a u ^ 7^q»2 + O ( ) 


■**Si ^^2 / P 2 “ “ ®2 ^2^2 “ ° 

(2.47b) 

k7*6- (Bj + K')7^(^-PC4 e " 0 

(2.470) 

and 


^2'^P2’S “ ° 

(2.48) 

and vector functions $^(x*t) and 

Batisfying 

the equations 


(uj. ♦ JD") 7^f + ( W 3 - 


- {a + & curl) ( ti - ♦ 2 * “ ^l’*! “ ° 

(2.49a) 

( U4 - l»") + « U2 + &0 "> ♦2 


‘^(oi curl) • $ 2 ^ **^2*2 " ^ 

(2.49b) 


If the functions are sufficiently 



( 2 ^ 50 ) 




differentiable# then w(x#t)# given by 

isfagrad <pj^ + curl v«grad tp^^fcurl 

p-(x#t)# and e(x#t) constitute a solution of 

4m 

(2*46) « Now the converse questions if the functions 
p^Cx^t)# w(x#t)# v(x.t), e(x#t) satisfy equations (2*46), 
are there scalar and vector functions cp^(x#t) satisfying 
equations (2.47, 2.48, 2.49) such that (2.50) hold? 

The answer is given affirmatively, and the representation 
(2.50) is complete. These results are given in the follow* 
ing theorem; 

Theorem (Atkin) 

Let p,(x,t) and Q(x,t) be scalar functions which are 

twice continuously differentiable on B and let w(x,t) 

and v(x,t) be vector functions whose fourth and third 

partial derivatives respectively are Holder continuous on 

B, th^ four functions together satisfying equations (2.46) 

in B. Then there exist scalar aftd vector functions 

9 . (x,t), 1 ). (x,t) which satisfy equations (2.47),(2.48) , 
i r 

( 2 . 49 ) such that w(x,t) and v(x,t) admit the represen- 
tations (2.50) in B. Moreover, it is possible to choose 
the vector functions ij£(x,t) so that the dilatational 
conditions 


div ^ 3 ^ ^ 
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are satisfied in B* 

With the aid o£ this theorem, the propagation of 
small amplitude waves in a non-»heat-conducting mixture of 
an isotropic elastic solid and an invidcid fluid was 
studied by Atkin [19]« Equating to sero the viscosity 

■ i 

coefficients %y (j*i# 3* 4), D»* and the thermal con-^ 
ductivity k# then differentiating each term of equations 
(2»47a) and (2«47b) with respect to t and eliminating 
and 0 by means of (2.47c) and (2.48), one obtains 

^ C |>2 - a(l-f) (*^ji - 4^) “ (2.51a) 

®4 ■^af (*i^^ ° *'‘^2 • (2.51b) 

^By the same process from equations (2.49a) and (2.49b), one 
obtains 

- a(l-f) (2.51e) 

a£(i(j " $2> “ (2.SXd) 

where ^ ” iP^ fractional contribution of the 

solid constituent to the mass of the mixture and 

°1 “ /?T Cjj) / P J . 

®2 ” <• ©2 P * ®a^ ^ 

■» (Kg + BiBj/P f " (l^+BiB2/p To^)/P2 


( 2 . 51 .) 



Eqtuations (2.5Xa) ana govern the propagation of 

small amplitude v/aves of dllatational type# that is# 
motions of the mixture in which the vectors w and v 
are Irrotatlonal# and equations (2* 51c) and (2.51d) 
describe motions of rotational type in which these vectors 
are solenoidal. Thus by the theorem (Athin) all suffi- 
ciently regular motions of the mixture can be decomposed 
into dilatational and rotational components* From equations 
(2* 51a) and (2«51b) one may obtain 


((V 


Pi 




^ 2 „2 




«* o 


( 2 . 52 ) 


where 


1 1 

,r ^ ^ 2.2 . y »^ 2 ^ 2 . 2 . 2 

y/2 

^ 11 

' 1 , 2 . 2 ,, 2 2»2 ,2 2 , 2. 2 

'' P 2 ” + «2 - < ^ ) 

1 

Vp 3 - (f(c^ e^) •> (1-f) (c| + c^))^ 

It has been shd^m that if Vp^^ and Vp^ real^ Vp^ 

is also real and the three dilatational wave speeds satisfy 
the inequalities 


Vp2 ^ Vp3 ^ Vpj_ 


(2.53) 
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The forsn oi equation C2«S2) suggests that at high frequen* 
cies there are snt^dlce dilatational wave propagation 
. associated with tb© speeds while at Iot fre- 

quencies, there is only one mode of wave propagation# 
associated with the epood ^h© oecond mode being a 

diffused disturbance « 

From equations (2.51c) and (2 .Sid) one may obtain 

2 2 

"“%■)+ a «0 (2.54) 

> 

The form of equation (2.54) suggests that rotational dis- 
turbances of the mixture comprise a single mode which has 
a wave- like character at all frequencies# the speed of 
propagation being of f^V in the limit when frequencies 
approach to sero and V in the limit when frequencies 
approach to infinity* 

So far we have reviewed the recent developments of 

the Interacting continue of an isotropic elastic solid and 
\ 

a viscous fluid* Due to the complication o^ the system# 
comparatively little progress has so far been made concern- 
ing the application of the linearized theory to particular 
physical situations# or the properties and understandings 
of the character of the system of the partial differential 
' equations • 



SufST^iary of thn Equations and Other Fonnalatlone 
A* Fully Goupls'l Mixture Theory 
At this point v.*e eurranarise tho pertinent equations 
«j^hlch oovsrn motion of a thermally and mechanically 
inteiTactinq continuous mixture according to the theory ex«» 
pour.if.ed l>y Atkin, Chadvick and Steel [18]. ife call this 
prob?i.om ^?y the name ” fully coupled mixture theory * ** 

8 train-d is placement equations 



(2.2S) 

rate of defornation-velocity equations 


*■ “<i. j)' “ ''(i. j) 

(2.12) 

vorticity-veloeity equations 


*ij “ ''[i. j] 

(2.27) 

‘"'h 

equations of motion 


* «j + Pi '■j “ Pi -at ' 

(2.29) 

— — . 

’"ij.i + “j * “ ^2 “St . 

(2.30) 


continuity equations 

- ^^2 - 



energy equations 


7 dt 9 Kan .10 mm ^ 


+ |’’'mm-«-|<Vm-Vm> /“O 


(2.31) 


constitutive equations 

®(ii) " ®1 " ®kk ■*■ ®8^ ^^2 “ 


P 


-5- 2(03^ -«- 2 u 3 * <2.35) 

’^(ij) “ ^*^2®2 ®2 "®8>®kk"P2“lO<®"^ 

P+P 5 - 

- < _J Oj + P2O5) (Pj Pg) ^4^kk ’** '^kk^ j 
P 

^ # <2.36) 

“i’- -^®l®kk.i* J®2 P2.i+“<«i"V <*•”> 

«i)pq<r‘pq"W 

®(ij] “ -’^Cij) ^ijp%“V"®"*^ij “ 

The complete initial*boundary value problem is sped* 
fled by the above equations ands the initial conditions 

*Obtalned from <2.31) by using (2*32) to <2.38) in <2.31)8 
Usually the form bo obtained is called heat equation. 
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(2.4I)« the boundary conations (2«42)# (2*4S)» and the 
material inequalities (2.39), (2*40) • The problem is 
solved if one can obtain at each place and t > O 

the functions and T* 

B. The Mixture Theory of Green and Steel* 

In a series of papers by Green and Steel [20], Green 
and Naghdi [14], Steel [21], a more tractable initial- 
boundary value problem than the fully coupled mixture 
theory of section 2.3 A has.been presented* The major 
difference between the linearized version of the theory 
presented in [15] and that of A lies in the constitutive 
relations. Green and steel’s relations follow from A if 
one sets equal to zero 

^1* ^3^ ^4* ^1^ Hi* (2«55) 

a", D, D”, 

C, Uncoupled Theory* 

By the term uncoupled theory we shall mean the 
initial-boundary value problem as specified in section 
2*3A\using the constitutive equations of Green and Steel 
presented in section 2«3B and in addition, neglecting 
the time rate of change of the dilatational effects of 
the solid and fluid components in the energy equation* 

If this is done, then the heat equation Is uncoupled from 
the equations of motion and the temperature may be treated 
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as a )aiown fiinctlon of space and time* 

2 •4. Single Constituent Theories* 

So that we can j^ke a comparison of the mixture 
theories presented in the previous sections with the clas- 
sical theories of elasticity and viscous fluids we record 
here the changes that must be effected* Our purpose is 
twofold: it allows us to drav/ upon the vast literature 

available in the classical theories of single constituent 
continuous media and it allows us to give meaningful in- 
trepretation to the mechanical and thermal material 
properties used in the mixture theory of subsections 2*3* 

Linear Thermoelasticity 

If the fluid component is not present then the body 
may be interpreted as a linear elastic solid undergoing 
thermal deformation in which the variation in temperature 
is small {22], (25], [24] . One seeks to obtain the compo- 
nents of displacement w, and the temperature T which 

A 

satisfy equations (2.25), <2*26), (2*29), in the absence 
of the fluid component and#with certain modifications of 
the heat and constitutive equations# satisfy (2*31) and 
(2*35). These mDdifications are formally equivalent to 
employing (2*55) and setting 

fiA# ^ (2*56) 


equal to zero 
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Equation {2.36) then interpreted as the cXassieal 
linear elastic stresses train law if we identify 

where Xg are the Lame elastic constants# y is the 

coefficient of linear thermal expansion and Kg is the 
isothermal bulk modulus# “ 3 * Prom (2*31) 

the modified heat equation becomes 

^^»n>m”PCa It" YTIC, “ O (2.58) 

mm ' e ot 13 ot mtn ' 

where 'p is the total density of the body# c the 
specific heat ^at constant a train « These coefficients are 
related to by 

a (2.59) 

• ' T 

The material inequalities (2 •39) and (2.40) simplify to 
^>0, c^>Q, 2yg+3Xg«=3Kg > 0. (2.60) 

A properly posed initial-boundary value problem of 

elasticity consists of finding w.‘ and T of class C? which 

% 

satisfy the modified equations and the following initial 


and boundary datai 
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® T a T T on B at t«0# 

v/^ sa w^ on for t > 0, (2*61) 

T a T “f © on ^ 

.where u^# # W^,©,P, and P^ are prescribed 

functions on the appropriate domains* 

We close this subsection with the remark that if the 

term 

® ^ ®nTjn <2.62) 

is ignored in '{2.58), then the resulting thermoelastic 
theory is known as the classical uncoupled thermoelastic 
theory. [22], [23] 

B. Linear viscous Fluid 
* 

If the solid component is not present then the body 
may be interpreted as a fluid undergoing thermal deforma** 
tidns in v/hich the variation in temperature is small. One 
seeks to obtain the components of velocity v^^ and the 

0 

temperature T which satisfy equations (2*12), (2*27), 
(2*30), (2*26), in the absence of the solid component, and 
with certain modifications of the heat and constitutive 
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equations <2*31)# (2.36}« These modifications are form* 
ally equivalent to employing (2. 55} and setting 

o# ' ^4^ ^5* ^3^ (2*63} 

equal to zero* 

Equation (2.37) Is then concerned with the unsteady 

linearized compressible flow about a state of rest of a 

heat-conducting viscous fluid# if we identify 
—”—*2 2 *— 

X# P ^ ^ 10 which X and ^ 

P 

are the coefficients of viscosity of the fluid# 
is the pressure of the fluid in the rest state# 

P (z: ^2 isothermal bulk modulus and p 

^is the product*' of the volume coefficient of thermal ex- 
pansion and the isothermal bulk modulus* 

From (2*31) the irodified heat equation becomes 

P ^ ■= ® ^2.64) 

where is the specific heat at constant volume# Y is 

the volume coefficient of thermal expansion, is the 

isothermal bulk modulus# and k is the thermal conductivity. 
.The coefficient is related a- by 

V / 

C a - 2 a . (2.65) 

P . 

The material inequalities {2.39} and (2.40) simplify to 
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3X + 2u > U > O, k > 0, ^ i ^ O 

> O and C^ i O. 

A properly posed initiaX«boundary value problem is then 
to find « p of class and T of class 

which satisfy the modified equations and the following 
init'^al and boundary dates 

A — A — A 

^i “ ^i* P«P+P# T»T + T onBat t»0 

V"p • 

on for t > o (2 *66) 

V « ^ 

T'o T + 0 on 

^Ic'k * ^ ^ i ^ 

where v^» ^ V^* 6# F 6^^ and. r are prescribed 

functions op the appropriate domains* 



CHAPTER III. RECiPROCXOT THEORS4S 

3*1. Introduction 

Prior to considering the reciprocity relation in the 
mixture theory# we will review th© well known reciprocity 
theorem of eXesticity. 

Suppose that an elastic body is subjected to two 
systems of body and surface forces. The work that would 
be done by the first system*© body and surface forces in 
acting through the displacements due to the second system's 
forces is equal to the work that would be done by the 
second system’s body and surface forces in acting through 
the displacements due to the first system of forces. Mathe- 
matically this is incorporated in the Bettl«*Rayleigh re- 
'ciprocal theorem. 

Theorem (Betti-Rayleigh) [25] 

Consider two equilibrium states of an elastic bodys one 
with displacements u^ due to the body forces and 

surface forces # and the other with displacements duo 

to body forces and surface forces • Then it follows 

that 

T.uf ds + j* dv ■ J* ds + J* F^Uj, dv . 

B ® 

A generalization of the reciprocity relation to dynamic 
problems is given as follows. 

0 

Theorem (Pung) [26] 

Consider two problems where the applied body force and 
che surface tractions and displacements are specified 
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differently* Let the variables involved in these two 
problems be distinguished by superscripts in parentheses 
such that the body force is the specified 

surface tract^.on is f|^^(Xj^;t) on , and the specie 
fled displacement is 9^^ (Xj^ft) on Where j • 1,2# 

Assuming that the action starts at t > 0 in each case, 
we have 

t 

J* J* (X,t-y)u^*^ (x.y) dy dv 

BO 

t 

+ J* J* (x.t-y)u|^* (x,y) dy da 

o 

t 

+ / J* (x,y)n^ dy da 

\ SB^° 

( 

t ' 

“ J J* (x#t-y)u|^^ (x,y) dy dv 

B 0 

“t 

+ J* / fr <x»t-y)uP' (x.y) dy da 

+ / J' (x,t-y)g|^^ {x,y)iij dy da 

3»B ^ 

As an Illustration of this theorem consider the follow** 
ing problems. By problem 1 let us mean the displacement 
and stress field in an infinite region that results due to 
the body force system 

o »(t) 



( 1 ) 

where is a fixed point in the medium and ' 

refers to a force magnitude. For problem 2 let us find, 
the displacement and stress field in the infinite region 
due to the traveling impulsive force system 

uc* 

xP O p<2) ~ 6(x2> SCxj) . 


!Chen substituting into the reciprocal theorem and using 
the properties of the Dirac delta function we find that 




(3.1) 


Pp^nj^<2cJ6(x3) dxj^ dx^ dx3| 6(y- ~)u t-y)dy 
x,/u 


V i 


* n ) *^1. ^ 

dxj^ . 


From this relation (p^pt) can be found when 
u. Mx, #0,0,t - is known. When F,' ' = 1 and 

= 0, Payton [27] has determined uf^^ by apply- 
ing th^ Laplace and Hankel transforms to the elastic equation 
of motion. The solution is 


X t 

ui^'(X,t) H ~ 0(x,,R,t), ui^^(x.t) » -i— P(R.t) 

^ 4ttR^ *■ 4itr^ 

^ (3.2) 

X.* X- t 

« - "■ ■V" 

^ 47ra^ 


where 
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P(R,t) =■ ^ H(t-R/Cj^) + ^6(t-R/Cj^) H(t-R/Cj) -^ftJt-R/Oj) 


a(x,.R,t) = - 1)“ H(t-R/Cj^) + “^Sit-R/Cj^) 


- ( 


^ - 1)| H(t-R/c^) ~ - 1)^ 6 {t-R/C 2 > 


1 

R a -5- x| •{■ X 3 ) ^ ; Cj and ©2 speeds of the 

propagation of dilatational and equivoluminal waves. 

Payton used the equations (3,1) and (3.2) to compute 


the u 


( 2 ) 


3.2. Reciprocal Relations for Mechanically and ther- 
mally Interacting Mixture 

t^e will investigate a reciprocity relation for the 
interacting continue of an elastic solid and linear viscous 
fluid using the theory derived in section 2 . 

We put 

^i “ ^i " ^2 “ ^2 ^ ' T « T + e (3.3) 


where is a reference position at time t » O, 

is a position at time t , P 2 is the density of the fluid 
component at T is the temperature at (x^,t), T 

is the initial temperature; then from (2.26) 




(3.4) 


all quantities now being regarded as functions of 
and t. Since initially- the medium is in equilibrium 
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under aero total applied force it follows that 


“i ** 


(3. 5) 


Let us consider the problems in which the body forces 
F^(Xi*t), , the specified surface tractions 

^i ' specified velocities are given 

functions of time and space# respectively# for solid and 
fluid which starts its action at t > 0 # with the initial 
conditions 


ww • 

v^eo#Ti«o, e-0 for tio# (3.6) 

Let the Laplace transform of a function u(Xj^#t) be 
written as u( 2 {^«P) where 

“J u(Xj^,t)dt . 

O 


^pply the Laplace transform with respect to the time t 
to every dependent variable. From (2.29)# (2.30)# (2.35)# 
(2.36)# (2.37)# (2.55)# (3.4) and (3.6) we obtain 


o , 
px#p 

"i + Px*'i 


(3.7j|) 



MM 


P*#P 

•"i Pz°i 

« PjPV^ 

(3.7b) 




Pl^2 ^ , 
P 



-{• 


+ a 


9 



(3.7c> 

(«4 - •|r^)%n®iH+2(o^+a5)aj^j^+(Og+ 

P 


(3.7d) 
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' To complete the problem let mb consider the boundary 


3B 

as the 

sum of the disjoint 

sets 

®®i 

and or as 

the 

sum of 

the sets (disjoint) 


and 

We specify 

for 

t i 0 

§ 

- V. - 

on 

3Bi 

(3.8a) 



(®ij + ’Tij)n^ - 

on 

3Bi 

(3.8b) 



«i . u. ,v^ » 

on 

®®1 

(3.8c) 



^ » T + © 

on 


(3.8d) 



q n « P 

^p p 

on 

®®2 

(3.8e) 


To aid in the computations we Introduce the following 
combinations of material constants: 


a 

P P 

- P 2 P ' 

^3 a Oi + t ^1 “ — ®2 " (3.9b) 

P 

p.ct« 

Yn = ?2 (-aft + -~) (3.9c) 

p 

Now consider two problems specified by equations (2.26) « 
(2.29), (2.30), (2.31), (2.35), (2.36), (2.37), (3.6) and 

(3.6) . We identify one problem by a superscript one on all 

field variables and a second problem by a superscript two on 
all field variables. For example # vj^^^ , 6^^^ will 

'satisfy (3.6) , (3.7) and (3.8) for , V ^ , 

q( 1)^ p(l)^ while solutions of (3.6), 

(3.7) , (3.8) ciifferent r|^^, SP* , uP^, Vp\ 

I 1» 2» 






. To derive the reciprocal theorem we begin witl> tjne 

equations of motion and the constitutive equations to which 

we have applied the Laplace transform and used the initial 

conditions as specified in equation (3.6). For the sake of 
« 

convenience these are 


O^V 

Pi#p 

- + 




(3.10a) 

--(j) 

pl#p 





(3.10b) 

^(j) 

gj ; ' S3 
X 

"ir 11 

p 

[<3) _ Ea „ S-tj) * 

. i ^ *’iwn, i 

Ct(Pw|^^ - V 

(j), 

i ’ 

(3.10c) 

^ ( T ) 

o :r ' =s 
ik 

^1 

P ^ik 

+ 2p S 

.2 snm xk 3 

«> . 


•« 

^(j) 

T . r 5* 

2.k 


# 



(B.iod) 

r ^2 
t- — a 

i ) T i 

2 “ + Yja”' 


®ik 

’^^rr ®ik 

■4* 


• 



(3.10e) 


for j 55 1,2. 

Multiplvir^ oquaicion i3,loa) for j = 1 by and 

X 

again for j = 2 by . subtracting these two results and 

then integrating over the region B, we obtain. 


J 

3 

*^( 2) 
u . 

X 

dv 

P^ * P 

r **^(2) 

Pi 

F^^^dv 


- ; 


^(1) ^ 




. (3.11) 

B 








--(1) 
u; ' 

^( 2 ) 

a'f' dv + 
P^#P 


Pi 

F^^^dv 


n 


B 




- 1 

* *•’ \ 
^i 

v> ' dv 



« 



B 
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\ 9 h^xe dv is the element of volume of B. 

Consider the first integral on the left and let us 
apply the divergence theorem to it* in this way we obtain 


rr(2) 




B 


^U) 

a'T' 

P>-#P 


_ r ^(1) j I* ^(2) '^(1) 

av - j » Si V® ” I “Ip Si 


(3.12) 


B 


where ds is the element of surface area of We note 

that the first integral on the ri t of (3*11) can be manlp** 
ulated in the same way and the re^ ^It is the same as (3.12) 
with the superscripts interchanged. 

Consider now the third Integral on the left* From the 
constitutive equation (3*10o) we may write 


J* u^^^ dv^’ra J vf^^ldv 

fe 'B 

as ■; P . P 

- j- pi;;! t - ^ s"’* 'i 5!" V 


(3.13) 


A similar result is obtained after the superscripts are 
exchanged* Using the constitutive equation (3*l0d)« we 
have for the second Integral on the right of equation (3*12) 




J i.P pi '’g 1,P pi 




J Pe' a^t'dv 


B 


ip pi. 


f . -I- 2^,S<^■hdv 

J ip ’■*^2 mm pi 


3^pi 
tfU) 


*• J- «pi 
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A similar result is obtc^ined for the middle Integral on the 
right of equation (3.12) after the superscripts are Inter- 
changed. When the equations (3.12) « (3.13), (3*14) along 
vrith the equations after the superscripts are Interchanged 
are substituted into equation (3.11) « we have 

J* + j* V® ^ av <3.15) 

B 3B ^ ^ B 

- J' f- ^ ^ ]n,d8 + / Pe^> 




P 


P 

( 1 ) 


V 

^(1) 


'mm -"2 
B P 


- X -F ®pi ^ M'"' V ■*■ 

a ]■ + J uf^'a^^^n^ds + J' Pav?^^*v^^^dv 

B , aa , " B 

r «r-(l)r ^2 ~<2) . ^1 '^(2) -i ^ r ^'^(l) ^1 ‘^(2)^ 

3B P B P 




( 2 ) 


B 

By they same process used to derive (3.15) from (3.10a)# we 
find that from (3.10b) 


r v^^^p^G^^^dv r ds J* Pa v^^^wf^^dv 

1 ^2 3. 1 P3L p r i 

as p ,3 * s p * ^ 

- XgV^fpt (- -T «2 •'■ 
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a. r r fl „ ~(2) . „ ~(2) . ,• ~(X) £2 „ ~(2)^„ 

^ j V + “ Cirri' jn.as + J v. . ^ ot-i®™ «v 

oB ^ —1mm p 2 i B P 


• r “■— G + Y P G 5 . Idv 

B 


ar(2>. 


(2) 


(3.16) 


Before we add (3.15) emd (3.16). we notice that the following 
volutae integrals 

p ~{2) ^2 ~(1) . ^ p „~(1) ^1®2 ~(2)j 

4 ''i. i J «l®m ^ { P’1 — 

- I' (- -~ a~ + p".iO, -0 ) 6 .dv 

J i,p ' P 2 *2 mm ^2 10 ^ pi 


B 


a. 


- r Pei^^[ —6 . + b. + a„e'-^'S„Jdv 

^ P P3- 1 pi S pi-* 

B 


(X) 


*.iay be recast by.Aasing (3.7) an<3 (3 *9) as 

\ 


B 


a, 

P P 


rp;?(2)r,^ , x^'yg) 

J Pn [ (-Og + _ _ ) p -JQ 


®~ ~(1) , 

- a,„a' ']dv 


(3.17) 


^ t i-s - ^ - v"' !■>'' • 

B P 

With the aid of (3. 17) the sum of (3«15) and (3«X6) can 
be written as 


(3*18) 


J + J' v?^^F^^*npde + j' “pj* »p< 3 s 

♦ L<»P’- 

• oB p p 

^ I < ^ <P2“l0^rr’- 


SB 
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“ .4 * 4 * 4b V" 

+ 4^uf>a^2)„^^3 

+ ; (yW - H- l^a2?[<2))nj,as 

3B P. P 

‘■^ / < ^ - «l^i^’>*‘m'*« + 4 ‘PaWri’ - 

JB ® 


Let the flux of heat and heat supply function be incorporated 
into (3 •18). The energy equation (2.31) reduces to 


d8 


+ a^d. 


■7 at ^ '*9-pp ■ P2 “ioSp ■ yS.P TS 


+ £ r 


O. (3.1d) 


Application of ^Laplace transform to (3.19) leads to 

P“7® ' PaWpp " I S.p ~ ‘ “ ° • <^-20> 

Consider the last integral on the left in (3.18). Using 
« 

(3.20) it becomes 


J- ■§r<l> (Pa,?<2). _ i + £ ?<2) 

B T T 

and if we now apply the divergence theorem to the qp p 6' ' 
term we may write (using (2.34) and (2.39)) 

K 


; Pa,et^>6<2)dv . I J 9<l)?<2)dv - I ; ds 


B 


T B 
' - ~i 

r B 'P ' P 


_ i j. ^(l)e<2)dv - ^ (u<2> - v<2))dv . 

7 ^ ''g »P ~ .n n n 


T as 

rC 

P 


(3.21) 
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We. note that a eimilar escpression can be obtained for the 
last integral o£ (3*18) by interchanging Indices* 

Sxibatituting (3*21) into (3*18) and employing the 
crane fomed boundary conditions (3*8) leads to the general 
reciprocal theorem in the Laplace transformed state* Be* 
fore giving this expression we introduce one additional 
condition* We set 

S<3> . (3.22a^ 

and require 

O (3.22b) 

P3. p 31 ' 

(3-220). 

on the boundary f dB^* This introduction somewhat simplifies 

4 

the notation but it must be recognised that only the total 
stress vector la specified on * 

l*hu8« by (3*8) « (3.18)* (3*21) and (3*22) we have 


/ »dv * f sr<«p,F<i>dv * f * r • 

B B SBjl 


- / vp>;?<l>n^a3 . r 


■+ J' (- ^ <*i( ^ ■*■ Px®2*P2^*“ 




J“ (v|2> - ^ Oj^( ^ - pj^h/pj^ * Px®2*P2^“ 

iff. I* 
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®®l 

^ p J.ir(l)g{ 2 )^^. 1 I ©<^>q< 2 S^. 3 ^ 1 j- 9 W?« 2)^3 

T B ? as, ‘ ^ SB „ 


' ~( 1 ),~( 2 ) „ ~( 2 )flv 


- ^ 7 1 ‘ 4 ’ (» 


T B 


»P ' P 


P 


« I + j‘ + J" Vj[^'gP*d8+ J* u?^^f^^ds 

8 B asj^ aBj^ 

»®i _ . 

+ J- (-R^^’H- ^ - p1^7/Pi + Pia2<P2^’ * 


dB. 


dB. 


Pn P 
P 


^ J (vf> » pJJ(l))[. . ^(2))/p^ ^ ^a,(p<2). i^)/^]n.ds 


+ J‘ + X - °l«m^7n„^d8 

®® 5 > 




dB. 


r 6 < 2 )?<l)dv-ij- 9 < 2 ) 5 ^ 1 )nds-iX ?< 2 )~( 1)^3 
T B as, as. 


. K' r ?( 2 ) J~{ 1 ) _ ~(l)ja^ 
' T '’b 'P P P 


<3.23) 


- Since the inverse transform of the product of two functions 
is the convolution of the inverse we obtains 



Let B be a bounded regular region o£ three-dimensional 
2?uclidean space occupied hy a mixture of elastic solid and a 
Nei^tonian viscous fluid undergoing a disturbance of small 
amplitude during the time interval t ^ O. We denote by ZB 
the boundary and by D the interior of B and introduce 
regions 6 of space- time defined by 

B ^ i (Prt) : p € B, t ^ 0), ^ =« { <p#t) * p € , t ^ O) • 


Wg use notation ZB^, ZB 2 and ZB^ , for arbitrary sub- 

sets of Zb and their complements with respect to ZB and n 
refers to the unit outv;ard normal vector field on 3 b. 

Sup^se .that tlio constants k, {x, a, k and satisfy 

3X + 2|i ^ 0 , 4 0 , e ^ 0, k i 0, k'^ i 4TOK , and that 

( % 
a-, 0-, ci^, Og, a,, and otg satisfy the inequalities 


>1 + “5 ^ 0 


+ y O5 ^ 0 


©2 


+ ag^O, 0,^0, Oj^( -=T - f-) 04 

p - 2p ^ 

“ «2 ^ ^ *°1 * “3* * “4 1 ®S* 


< = “2 ^ V • 


Let the mixture of elastic solid and viscous fluid be sub- 
jected to two systems which are distinguished by superscripts 
in parentheses. Let furcticns 

class and of class on B which satisfy 

equations (2.26), (2.29), (2.30), (2.31) on B* , equations 
(2.25), (2.27) and ( 2 . 34 ) to (2.33) on B and the 
subsidiary conditions 


62 


^ • Or ■ 0, « O, ti w 0, 0 » O on B at t • 0» 

i • ®pi "p ' ’^pi "p ■ «t » 


-iJS - sP’* 


on 0B. for t ^ 0 


U. « o|^' ♦ V. a vP* on 0Bj^ for t2iO, 

y(j) « T + on SBj, q^^^n a on SBj for t ^O, 

where ,f|^^ ,0^' ,P^^' .P^* .G^'and r^^* 


are prescribed functions on the appropriate domains and 
p^r ^ 2 » ^ given, positive, constants* Then the work that 
would be done by the first system in acting through the 
velocities of the second system and the work that would be 
done by the second system in acting through the velooitiee of 
th^ first system satisfy 


t • t 

?2 / J <x,t-y)oP’ (x.y)dy dv+p^J J uP' (x.t-y)pP' (x.y)dydv 

B O BO 


t t 

*i S vPNx,t-y)gPNx,y)dyds+ j* J «PNx,t-y) fP^ (x,y)dyd* 
3b/ . . te^O 


t u 

+ / J vPNx,t-y)ir^P (x,y)npdydSf / J oP' (x,t-y)ff^P (x,y)np4yd» 

8B, ° as- ° 

t **P 

- j*\ J rPNx,W)[ ~ o^Cpj-pp' (X,y)) 


* P,a 


PtP 


+ (x,y)-p2)l <3y + 

P 


^Attention is called to the fact that only S. is known on IM, 
and hence only the sum + 9^, is given* 
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* L Ja (x.t-y)) [- ^ o, (p,-p<^J (x,y)> 

sa, pjp * * * 


+ <pj^’ (x.r) - Pg) ay 

P 



+ £ J* X <x,t-y) (x.y)dy dv - 

T B 0 

' I Jb { [q^ 2 ’^x,t-y) ©<i) (X,y) }np dy ds 
\ 2 

* I '!' *^0 <x.y)dy ds - s' J- Hu< 2> (*,t-y) 

* 2b, T b 0 P 

- (x,t-y) ]8 ^i^ (x.y'dy dw 

1: 

“Pa XX v^^Nx.t-y)G? 2 > (x.y)dy dv 

B 0 

^ t 

Pi J X «i MX,t-y)Pp) (x,y)dy dv 

‘b 0 

t 

+ X X Vi Mx,t-y)gp) (x.y)dy ds 

SB ° 

1 

t 

+ X X «1 <»*t-y)f< 2 > (X,y)ey ds 

2B.0 

X 
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+ J / (X,y)n dy ds 

OB- 


+ X ,)■ «i^^<x,t-y) (x.y)n dy da 

8B, ^ . 


J’ j’" (x,t-y[~ a^(pj_ - p<2) (j,,y)) 


^'1 ° 


PlP 


/ 1 /3\ _ 

+ r- Oj (Pj <x*y) - p)] dy da 


+ / J* ( (X, t-y) -u (X, t-y) ) { - ^ o. {p. *p 

affj^ ® PjP 

_ ®2^p2^^ ” PaJl “i *^y 

^ * * 

- / J‘ n dy da 

;b, ° 


* L L £p2®2''m^’<’‘*<=-yj - «i<^'(x.t-y)] dy 


( 1 ) 


aa 


mm ^ 

+ £ J J (x,t-y)0^^^ (Xiy)dy dv 

T B 0 

" ^ J* / <3o^Nx,t-y)0<2^ (x,y) n_ dy da 
*' aBjO ^ ** 

• i J / P^^Nx,t-y)e^2) yj^y ^3 _ ^' . J.*" 


- <■ 


T B O 


- (x,t-y))6^^^ (x.y) dy dv 

F #P 


(x,y)} 

ds 

(«p^^ (x,t-y) 
(3.24) 
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. 3«3« Special Cases • 

A* ReciDroeity Relation for Heat Conducting Mixture 

*ii iiiiiin m iim n i — Ofccii— ^itgaamtatii— uti— Hw— — — — i— 

of Linear Ela stic Solid and Non-Newtonian Viscous Fluid 

Let B b3 a bounded regular region of three-dimensional 
Buclidean space occupied by a mixture of elastic solid and a 
non-Newtonian viscous fluid undergoing a disturbance of small 
amplitude during the time interval t ^ O. Yfhen the viscox^i 
fluid is non-Ne\'/tonian, the restriction (2.55) does not hol4 
in general* Let the mixture of elastic solid and non-Newtoaian. 
viscous fluid be subjected to t^o systems v/hich are distin«» 
guished by superscripts in parentheses, where each of these 
systems has same initial and boundary conditions stated in 
the previous thporem* Then the reciprocity relation is again 
an integral relation obtained by adding the follotying equality to 
(3,24). 

( 

- 2 J’ I'*' fiipV«x.t-y)(X3f^J^(x.y)6p. + 2 (i3f^^><x,y)) ay dv 

B ® 

+ a'J- vJ 2 )(x.t-y)€^pg(r^J>(x.y)- A^>(x,y)) dy dv 

B 

■ X X (x,fe„y) (X^a^i) (x.y) (x.y)) <3y av 

• B O 

- - a' J (x.t-y) (x.y) - (x,y)> dy d» 
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- 2'X + 2i*3fj?(x.y)) «!y dv 

B O 

♦ a' X (jc,t-y) 6jpq(r^** (x^y) - (x,y)) dy dv 

BO' 

B ^ 

< 3 . 25 ) 

B. Recipyocitv Relation for Mixture of Linear BXaatlc 
Solid and Newtonian viscous Fluid in Isothermal Process 

Let B be a bounded regular region of three-dimensional 
Euclidean space occupied by a non-heat conducting mixture of 
elastic solid and a Newtonian viscous fluid undergoing a 
disturbance of sjnt^all amplitude during the time interval t ^ 0* 

■ 4 

Let the mixture be subjected to two systems which are dis- 
tinguished by superscripts in parentheses, and let the process 
be isothermal* Moreover, except the temperature terms, we let 
each of these systems have same initial and boundary condi- 
tions stated in the previous theorem* Then the reciprocity 
relation is given as following; 


Pj r r* vp* (x,t-y)GP^ (x,y)dy dv 

BO 


' w 

**“ Pi J / ttp)(x,t-y) (x,y) dy dv 
B ^ 

a 

♦XX vPMx,t«y)gP' (x,y)dy do 




♦ X X «pNx.t-y) fi^Mx.y) dy do + 


Cl) 
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+ J'_ X” (x,y)i»p dy 




pi 


ds 


^ t \ 

+ X_ L »i ' (x.t-y) (x,y)np dy da 

- X X^ (*,t-y)t Oj^ip^ - p<l> (x,y)) 

fta- ^ P,P * ^ 


fts, 

« • ' 
x.a 


n;* /«(!) 


<P2 ' (x,y)-p2>]nj^ dy da 

•" 

♦ X X (X. t-y) -op> (X, t-y) ) t - ^ a, {p, 
33.0 P^P 


'X 
P»« 


+ “^ (p(l)(x,y) . Pj)]n^ dy da 

- X X ?2“2 *s5« «y a« + X X (Pa< 

' 33, O 

■ (x.t-y))nj„ dy da 

■ ?2 / <L (x,t»y)GP^ (x,y) dy dv 
B O 

•» ?x / (x.t-y) (x,y) dy dv 

♦ X X^ (X.t-y;g<2) (x,y)<Jy as +J> A d) 

8B, o 8B O 


♦ X^ X^ (x.t-y)v^f) (x,y)Bp dy da 

3B 

■•■ X X Ot^Nx,t-y) OoP (X,y)n_ dy da 
VS ^ ■ 


•Pj[^* (x,y> 


2 v^*>(x,t.y) 


(X, t»y) (x,y) dy ds 
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“ I J’ (x*t-y)C- Oi<?i * (*«y)) 

'®i° PiP ^ * 

« » ♦ 

•*• • ^ ®2*^2 ' - p2>ln^ dy ds 

P 

"t ’p 

■*■ (_ "f (vP' (x.t-y)- (x,t-y)) (- o,fpr.*pP^ (x.y)) 

dB O *■ ^ p,p ‘ ^ 

«» 

+ ^ ° 2 ^** 2 *^ “ P 2 )J“i <3y <Ss 

P 

* / J1 P2“2 (x.t-y) dy da ♦ | J X PjOj (x.t-y) 

as, ® as, 0 

1 1 

(*'**’5r)l n^n ^y ds 

( 3 . 26 ) 

C. Reciprocity Relation for Heat Conducting Mixture 
of Linear Elastic Solid and Newtonian Viscous Fluid Occupying 
Infinite Region 

Let three-dimensional Euclidean space be occupied by a 
mixture of elastic solid and a Newtonian viscous fluid under- 
going a disturbance of small amplitude during the time interval 
t ^ 0* Since there is no boundary in this infinite region 
casop Regularity conditions will be considered in the place 
of boundairy conditions. Let the mixture be subjected to two 
systems which are distinguished by superscripts in parentheses# 
where each of these systems has same initial conditions stated 
in the previous theorem and satisfies the following regularity 
conditlonss 
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9^^^ ■ ■» O as (icl • foe i •> 1,2. 

Then tho raelpeocity thooeem 0 , 24 ) eeducea toi 


hS 

B 

!' 

0 

(x#t*y)G^^^ (x#y) dy dv 

mm 

* Pi 

J* 

J uP^(x,t*y) (x,y) dy dv 


3 

0 

+ 1 

! 

D 

1 r^^^ (x,t*y)0^^^ (K»y)dy dv * 

T 

0 ■ 


T B 


O ‘ P 


) 0 ^^Mx,y) dy dv 

P #P 




h -f L <x,t-y)G^‘ 2 » (x.y) dy dv 
B 


( 2 ) 


^1 ^ J'o (x»y) < 3 y dv 

B 


+ £ J* I 0 ^^Nx#y)^y fiv - ~ J f [ uj^^^(x,t-y) 

? B o T B o P 


(x#t-y)l 6 ^^^ (K,y> dy dv 
P • *P 


( 3 . 27 ) 


where the integration on B is over the entire three* 
dimensional space, |x| < « # |y| < ^« |e| < 


D. R^i o r peitv, Relation for Heat Conducting Bias tie 

* SpUd 

Let the elastic solid be subjected to two systems which 
are distinguished by superscripts in parentheses. Let the 
functions and be of class on B and the 
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subsidiary coKiditions ar© 

O9 Q ® o on at t =* 0^ 

^Cj) „ ^ ,j, (9^3) Qjj ** 0 *^ ^^2 ^ 

bliQre fP^, are - 

proscribed furictlOs^s on the appropriate domains, and 
‘2’ are given positive constants. Ihen by the result 
ot Section Z*4 Part h, it follo-^s that the work that 
would be done by the first system in acting through the 
displacauiont of the second system and 'the work that would 
be done by the second system in acting through the die** 



placeoent of tl^a first syctcim satisfy the following re- 
la t^n. 


p| I (x.y} ay dv 

B O 


■*"11 ^ss,t-yj iK,v) dy ds 


d3 


O 

X % 


■*■11" (ss,t«y)a4^ (x,y)n dy ds 


Cl) 


53 


O 


P3. 


P 


1 


C2) 


( 1 ) 


< £ J r^‘^Mx,t»y)0'*' Cx,y) dy dv 
s o 

“ r j' I 0:,t-y)©^^* (x,y)n„ dy ds 




p 


" h s S «J5,t-y) 8^^^ (X,y> dy ds 


■i O 


4si 
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P' J J wi^Nx.t-y) (x,y) dy dv 
B o ^ *■ 


( 2 ) 


■•• X X w|^Nx»t-y) (x,y) dy ds 


( 2 ) 


dB. 


0 


+ X X (x.*-y)Opi' (x,y)np dy ds 
° 


( 2 ) 


♦ ^ XX* x*^Nx,t-y)»^^Nx,y) dy dv 
^ B o 

" — X X* (x,t-y)0^^* (x,y)n_ dy ds 

- 3B, O ** 

- i X X* (x,t-y)6**^ (x.y) dy da (3.28) 

as, o 

2 

We remark that this result agrees with the well known 
leclprocity relation for elastic solid in isothermal case^ 


£• Reciprocity Relation for Heat Conducting Viscous 



Xiet *the viscous fluid be subjected to two systems Which 
ire distinguished by superscripts in parentheses. Let the 

e 

Unctions vP^ be of class and be of class Oh 

$ and the subsidiary conditions are 
B n “ 0, 6 « O on B* at t » O # 

^pi"p “ sP’ o« '^i “ v«> on 8ff^ for t i 0 , 

a T + on SEj, ^p^^n “ on SBj for t i 0, 

where P*^^, and r^^^ 

are prescribed functions on the appropriate domains and p »T 
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are given* poeltive* constants. Then by the result of 
$ection 2«4 Part B* it follows that the work that would 
be done by the first system in acting through the 
velocities of the second system and the work that wuld 
be dene by the second system in acting through the veloc* 
it lea of the first system satisfy the following relation! 

? r (x.t-y)Op^ (x,y) dy dv 

BO 


+ J S (x.t-y) (gP' (x,y) + p dy ds 

M 1 

♦ / I vp' (x,t-y) (iTpP (x.yjHp ♦ p 02"i* 

1 

V 

•*•—//* (x,t-y) 6 ^^Mx»y) dy dv 

T B Q 


(X) 


* - J’ X (x,t-y)©*^> (x.y)n_ dy ds 
- ~ X X P^^V<x«t-y) 9^^^ (x,y) dy ds 



vP^ (x,t«y)6p^ (x,y) dy dv 

vp' (x,t«y) (gpMx.y) + "p By 


ds 
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<j5,t-y) («»y)np + p 

t 

j* (x#y) dy dv 


t 


X C«.t-y)©<^^ (K.y)n^ dy ds 


X (x#t-y) (x,y) dy dB 

0 


(3.29) 


;ion in Mixture 


Theory 

Suppose that infinite three-dimens lonaX Euclidean space 
is occupied hy^^e mixture of elastic solid and Newtonian 
viscous fluid undergoing an isothearmal disturbance of small 
amplitude during the time interval t ^ O* X*et the 'mixture 
be subjected to two systems which are distinguished by super- 
scripts in parentheses with the following specified body force 
systemss 

= aP6(pj^)6(t) 

*0 

where p. is a fixed point in the region and a^^' refers 
to a f^rce magnitude, . 


Where p^ 


p(2) ^ Q 

6{pj^) 5(t) 

is the same fixed point in the region and 
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refers to a force magnitude* hek us assume that the veloc* 
ity fields of these preliminary problems are known# that is^ 

known. Suppose that the 

mixture is subjected to an arbitrary body force system 

and . 

^en the velocity fields and at the point 

p^ at time t are given by the following integrals as U 
result of (3 .28): 


^1*1**^ (Pl.t) * Pj J' J" (x.t-y)oP^ (x,y) dy dv 

BO . 

+ ?1 / J" (x.t-y)pP^ (x,y) dy dv , 

a o 

\ 

pap>vP^p^,t) « Pj J_ vP>(x.t-y)Gp>(x,y) dy >dv 


t 


+ Pi S Jf (x.t-y)pp’ (x.y) dy dv . 


<3) 


S O 


(3.30) 



CHAPTER IV. 


. A FUNDAMENTAL ONE-DIMENSIONAL 
INITIAL-BOUNDARY VALUE PROBLEJl 

4.1. Introduction 

In linear therrooelasticity, for one dimensional model, 
a scries of papers by Danilovskaya [28], Sternberg and 
Chakravorty [29], and Muki and Breuer [30] have answered 
some of the basic questions concerning the .effects of the 
inertia terms in the elastic equations of motion and the 
effect of the mechanical coupling terra in the Fourier heat 
conduction equation. 

Recently, Martin [31] studied the initial -boundary value 
problem corresponding to Danilovskaya's in the framework of 
the linearized interacting mixture theory. In [31], [32] a 
mixture of linear elastic solid and viscous fluid occupying 
a half-space undergoing deformation due to a transient tem- 
perature on’ the boundary was considered. Method of the 
solution was that a parameter occuring in the diffusive re- 

y- 

sistance vector was- used as the basis for a perturbation 
procedure in the equations of motion. 

Here, we consider the same problem but use a different 
method of solution. Let the mixture occupy a half-space 
and let its motion be restricted to one space dimension. We 
‘prescribe a step function temperature on the face of the 
half -space whore the face is constrained rigidly against 
motion. 
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4«2» Formulation of the Problem 

Departing from the indicia 1 notation# we describe the 
cartesian coordinates as (x,y, 2 ) and consider the mixture 
of an elastic solid and a viscous fluid as occupying the 
region ^ >0 [31]* We assume that the mixture is subjected 
to a time-dependent temperature field of the form 

T * T(x#t) <4#1) 

and is restricted to uniaxial motion, so that the displace- 
ment vector of the elastic solid has components 


W. « (w(x#t)r 0# 0) (4.2) 

and the fluid velocity vector is 

V. = (v(x,t), 0# 0). ' (4.3) 

It is convenient to define a nondimens ional temperature 
field in place of T in (4.1) by 


0(x#t) as I ■ 

T 


(4.4) 


where T denotes the temperature of the equilibrium state. 
Substitution of (4.2) and (4.3) into equations (2.12)# 
(2.-2S) and (2.26) shows that we may write 


” dx ' ^x dx ' St “ dtox ' 


(4.5) 


2>w, 


p^(x.t) = p,(l ~ 


<4. 6) 



77 


Bp 5 

at" 


-J* p 


5v 
2 3x 


= 0, 


(4.7) 


as the only non-vanishing kinematic relations# and we note 
that all other strain# rate of deformation and vorticity 
components are identically equal to zero. The constitu- 
tive equations (2.35) to (2.37) under the restriction (2.55) 
become 


p 


-2. 

3 w 


^1^2 


Bp, 




(4.8) 




oW 


a^(x#t) « ( 2 (a^ta 5 ) - (— - a^) ) 

H 


a. 


\ 


+ O^T0 + (— 4- Og) T1 + # 


(4.9) 


ay(x,t) = a^(x#t) = 




3w 


a, 


“ q.)-^ + agT6+ (~+ag)ri + a. # 

p ^ ^ p ° ^ 


a *= a » a « o# 
xy xz yz 


(4.10) 

(4.11) 


ir^(x.t) = (2u+X)g+P2(^-ag)g-P2a^^T9 

% P 

P+P 5 - 

“ ( ' ^2 ^2^2 * (4.12) 


Ty(x,t) -T^(x.t) =X^+p2(^ “2“«8^-^-p2“lo’'® 
P + P 5 _ _ 

_ ( — _ — + P2“e^ " P2“2 ' 

P 


(4.13) 
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The equations of motion (2.29) « (2»30) and the energy equa* 
tiOA (2«31)i under the constraints (4*2) to (4*4)* become 


■sr “ «x 


air 

: 

Sx 


p ^ 




^ “ ^2 3t ' 


fl OtgT+K* -2 ?^v 

“7'^ + < — = — > a3r=°- 


(4.15a) 

(4.15b) 

(4.16) 


To complete our formulation of the initial«boundary value 
problem we prescribe that for t < O 

w(x,t) = |^(x,t) = v(x,t) = 0, 

P2 = Pjf S(x,t) * °* (4.17) 

In addition, we require that on the boundary x = O, 

%(b.t) = h(t)*, 1^(0, t) = 0, v(0,t) - O, (4.18) 

while as x-«®, we stipulate that 

0{X,t) ,W(x,t) , P 2 (X*t), V(x,t), 0^(x,t), 7T^(X/t)* Cy(x^t) 
and \iry{x*t) approach to zero* (4.19) 

At this point, we introduce dimensionless variables. For 
- 'this purpose we use the notation introduced in (3.9) 


4r 


h(t) is the Heaviside unit step- function defined to be 
zero for t < O and one for t >0* 
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wherein the relation (3.5) is used* A direct substitution 
of ( 3 . 9 ) into (4.9) yields 






(4*20) 


which^ if the material were elastic, would lead us to ex- 
pect ^ role of the Lame' constants 

(Xg + 2ug) while a^T would play the role of 
(2u„+ 3X«)ci„ where a„ is the coefficient of linear 

£1 £# £« £# 

thermal expansion of an elastic material. With this in 
mind we choose a velocity , 





^2 *** ^^3 


(4*21) 


V > . . 

which would be the irrotational velocity of sound if the 

material were elastic. Since < O by (2.40), we define 


2 



(4*22) 


By a dii^iensional analysis we have that c^ is a velocity 
while (4.22) has dimensions of length squared per unit 
time. Thus, if we take 



(4.23) 


then a dimensionless x-coordinate and time are given by 


r 21 - 
^ * a ^ 2 

U) 


C^t 


U) 


(4.24) 
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Proceeding further, wo introduce non-dimensional partial 
stresses, solid displacement, fluid velocity, densities 
and diffusive force by 


A ^ 

''x ~ + 2{3. 


A w 

# w = r* 


ct, . 0 

i— . $ i: 

0 ^2 + 2(53' ‘^y P2 + 2P3 


A 

=S 


A 

A X 

+2(3, 


^ V t^ 

V a 0 

a 


A 




X 


Pj + 2;i3 




A 

“ pj + 2P3 


(4.25) 


In addition, the following quantities are conveniently 
grouped : 


^2 2utx... a a 

tQ(P2+2P3> P2+2P3 ^2 ^^ + 2^3 

. £2 , "^2 A "^2 Pl °‘ 2 . 

■p ' ^1“P2-" 2P3 - ^'^0-p2+2P3 <^1- - > 




aaT+K» 



Pjai,,T.K- 




(4.26) 


Incorporating all of these changes leads us the following 
summary: 

Constitutive equations 

\ 

Ojj(C» t) = t)+[ (1-f) 0 ^ - 6j^]t!2 (C» t) 


(4.27a) 
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^ " “®o + fa^) 5 Y . (^j . TI . _ 


a^efC t) 


+ (63- (l-f)Oo>’T2<^' 


w ( C » f ) -~£^q 

SC 


&rj. 


at 


o , a*) 


♦ =r^ t|^(C.T) - V(C.T)J. 


(4.27b) 


(4.27c) 


"07^20, +dj3(CT) 


+[{X-f)o^- 6j^]il2(C. T) + a^. 


(4.27d) 


Ojs(C'T) s ay(C,T), 


3v - -A 




(4.27e) 




+(62- (l-f)®^,]T)2(C.‘r) -o^ , 


(4.27f) 


IT (0t).^= 


(4.27g) 


Equations of motion 
,2.a at^ 


2., 55,.^ 


It" ‘1 it 




2 ' 


-2, ^2A 

bT ^1 ac^T *^2 ac ' 


(4«28) 




(4.30) 



' Initial Conditions 






o. 


(4 >31) 


Pot t ^ O# we take 
A ^ A 

w(C,t) <«v{C, t) -TIjCC.t) • e(C,T) • 0. (4.31) 

Boundary conditions 
At ( m o# we take 

A . 

9(0,T)»h(T), -^0 ,t)«0, v(0,t) -0. (4.3:) 

As C •« a. we take 

0(C.T). w(C,T).^ v(C,t). T12(C.T),$^(C.T), 

^ J 

and (4 *14) 


For this initial**boundary value problem specified by equa«* 
tlons (4e27) to (4.34)/ we assume that the initial street 

of the solid constituent is zero# ise*# 

■% 

« 0/ (4.3S) 


that the mechanical coupling terms in the heat conduction 
equation (4.30) can be neglected 


6-1 « 0, e 


O, 


(4*36) 


and that the constants of the mixture satisfy 
9 >r »- ^2 » t >0, - &2“^1 ^ *2 


i 


(4.37) 
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Where t 



and 



(4.38) 


The assumption of (4,36) is analogous to neglecting the 
effect of the mechanical coupling term in the Fourier heat 
conduction equation in linear therraoc las tic ity theory [22]# 
[23] and is justified by the conditions (a) and (h) of 
Section 2.2. 

Various material constants v;hich appear in (4.37) have 
to be determined experimentally for the mixture. Ihe re- 
strictions on the material constants# (4^7)# are the results 
of (2.39)# (2.40)# and their interpretations in Section 2.4# 
Single Constituent Theory. For additional references consult 
[33] to [36]. 


With the aid of (4.35) to (4.38)# the equations (4.28) to 
(4.31) are written 

(4.39) 






2 ^ 


1 ac" 


dT)- 


2 A o^ 2A 3 jiA 

. 0 w ^ ow ^ 2 o V .A J . c z OV (A 

6, + t *5^ + s — ^ - tv - * 5 T? + 6^ - 57 T r ^ # (4.40} 




'2 2 ac 


8^0 








dC 


= o # 


=s o . 


(4.41) 


(4.42) 


4.3. Solution by Integral Transforms 
we denote the Laplace transform of a function P(C#‘t) 
with respect to t by F(C#p)f where 
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F(c#p) f dT 


(4.43) 


c+ic3 


2iri P(?,t) e : P(C,p)e^'’^ dp, 


(4.44) 


C-2L«> 


In (4*44) c is a positive real nuniber such that the path o£ 
integration is any vertical line to the right of all eingu* 
larities of F(C#p) p-plane. 

The solution to the thermal problem* equations (4«4X)* 

(4<32) to (4.34)* in the <C»P) pl^ne is 


F(C#p) « ~ exp(-p^5) 

P 


(4. 45) 


Novf>^trans forming the equations of motion (4.39), (4.40) and 
the equation of continuity (4.42)# then combining these 
equations with the solution to the thermal problem (4.45) 
gives 

6 d 

(D^ - pt - p^)w+ (“ n^^t)v»-^ expC-p^O (4.46) 

P . pB 

+pt)w+ [ (s^ - t-pr]v« y exp(-p^5) (4«47) 

P 

v;here the differential operator is defined by 




D « 


A. 


Due to equations (4.32) to (4«34)* the transformed boundary 
conditions and the regularity conditions are 


v(0,p) » O* v(0#p) », O 


(4*48) 
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and as ( 

w(5,p), v(C,p) -«0* (4.49) 

solving for the displacement of the solid constituent and 
the velocity of the fluid constituent from equations 
(4.46) and (4.47} for the homogeneous solutions only Which 
conform with (4.49). we have 

w^<C»P) “ B{p) exp<-|^ (Sl + STjlO 

r^* 

+ D(p) exp(-J~ ^^1*^2^^^^ (4.50) 

4 

v^ = [p{(t+p) (ps^-dj) +t 6 j) *•* 

.« 1 

+[p( (t+p) (ps^- 62 ) +t 6 j) -|p(gj-g 2 )^]D(p) 

exp(-||- (gi-gj)?)* (4.5X) 

1 

where B(p) and D(p) are integration constants# and 
9 ^(P) = (ps^-^j) (p+t)*h(pr+t)-«- 2 t 5 ^ 

0 

•f 2p^(ps^-6|-62) ^(py+t+tr) (4.52) 

92 (P) = (Ps^-62) (P+t) + (pr-}-t)+2t 

- 2p^(ps^-6^-62>^(pr+t+tr) (4.53) 

f|(p) =*ps^- * 1 “ ^2 • 


(4.54) 



ss ^ 


Hero we hava acsumed that 


>0 and Re(r^^‘ (gj-gj)O >• 6 (4.55) 

for p such that Ro{p) s- !l s- 0, and vfe will show this to 
be so later* Let us indicate particular solutions to the 
equations (4.46) arid (4.47) by a subscript "p“. Then we have 


Wp(C,p) " ~ expC-p’^* C) 
Vp(C.p) = 


(4.56) 


(■ 


p{6i+6a-ps^+t6i + (ps 2 - 62 ) (p+t) )~ 




p(ts*+r 6 i )+t( 6 i- 62 ) 


■ps^dy-iadi+dida 

— w^w m ■ 

1/ 2 

p(ts*+r6i)+t(6j-6, ) 


f 


.V 


2 0 


where 


(4.57) 


(4 .58 ) 


Ai = di(p(s 2 "r)- 63 -t)+d 2 ( 6 i+t). 

As ® PP^'2(ps2-6f"62-t"pr-(t+p)(psa-62)-3t6i+p2r+pt+ptr) 

(4.59) 

The entire solutions are 

^w(C.p) = Wjj(i;,p).^+ Wp(i;,p), (4.60) 

v(C,p) = Vj^(C.p) + Vp(C,p). (4.61) 

TO detertaina B(p) and D(p) wo substitute equations (4.60) 
and (4.61) into the boundary conditions (4.48). Then 


+ -vP; ■ A„^p) 


(4.62) 
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(gi+92 )* 

B(p)pt(p+t)(p8®-6g)+t6i- — — ] + 

{oi -g<i 

D(p)p( (p+t) (ps3-6a )+t6i- -) 

+ (p(ai+6s-psa+t6i+{psa-63)(p+t))^ 

"2 

(ps®-6a)di dida 

+ -■■ + “i77^ " ® * (4.63) 

Solving for the tinlcnowns b(p) and D(p) from equations 
(4.62) and (4.63) simultaneously, we have 


b(p) 


*1 +*2 -ps® + |( (pe3 -6a ) (p+t )+(pr+t )+2t6i )) 


{{ps®-5a)di + 6id. 


I I 

'pp"'* 


gi9s» 


^31 

2A2 


(4e64) 


b(p) ■> I6j+6a-ps*+ |{(ps8-6a)(p+t)+(pr+t)+2t6i}l 


{(ps2-6a)di + *1^2 I 

pp^'®gig2 


2A2 


(4.65) 


The Laplace transformed displacement of the solid constitu- 
ent and velocity of the fluid constituent are now explicitly 
given hy equations (4.50), (4.51), (4.56), (4.57), (4.60), 
(4,61), (4.64) and (4.65). These quantities constitute the 
complete solution, in the transform plane, to the initial-’ 

, boundary value problem posed by equations (4.32) to (4.34), 
(4.39) to (4.42). Stresses of each constituent may be ex- 
pressed inunediately in terms of the transformed displacement 
of the solid constituent and velocity of the fluid 
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constituent toy means of the constitutive equations (4«27a) 
to (4,27g)* Now that the Laplace transformed displacements 
of the solid constituent and velocities of the fluid con* 
stituent are in their simplest form \te proceed to invert 
these expressions • 

4«4« Inversion 

A* Location of Zeros of gx(p)g|(p) 

As a first step toward the inversion of w((;,p), w© 
examine the multiple valued functions appearing in equations 
(4*52} and (4*53)* Set 

- 63 

^ 

ea = . (4.67) 

Then equations (4*52) and (4*53) become 

gi(p)={ps®*^a)(p+t)+(pr+t)+2t6ji+2r^/2sp^/a(p+€i)^/a(p+eg 

(4.68) 

9a(p)*‘{ps*"*a ){P+t)+(pr+t)+2t6i-2r^^*sp^'*(p+€i )^^*(p+ea )'^® • 

(4.69) 

We teke the domain of definition of p^^^ as the entire 
p-plane cut along the negative real axis, and we choose a 

' 1 /« 

''branch of p through the requirement that 

p^/a ss ^ for p “ jg > 0.*^ 

*</T refers to the positive root for real positive £* 


(4*70) 
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1 / 

Wo take the domain of definition of (p+€i) as the entire 
P'^pXane cut from ->€x to -OD along the negative real axia« 
and we choose a branch of (p-fci) through the requirement 

;; 'W.: T;. . . 

that 

(p+Cj.)^/® = for , p = f > -€x ** (4.71) 

Wo take the domain of definition of (p+€a)^^^ as the 
entire p-plane cut from -<-2 to --oo along the negtive real 
axis^ and we choose a branch of (p+G 2 )*^^^ through the re- 
quirement that 

« .^7+^2 for p =-* i > -Cg (4*72) 

tJow we will choose a domain in which 9x(p) ^nd ggCp) ere 
singl^Walued . In view of (4.70), (4.71) and (4.72), we 
find that 9i(p) does not have branch points but ggCp) does* 
To see this we note th^t all possible branch points of 

9x(p) 9a (p) zeros of gi(p) 92 ’(p) which is a 


fourth degree polynomial, 

2 (ts*-r -62 ) 


gi(p)92rp^ “ + p®(- 




)+p®t 


2t(-l+26i-6j|-2r ) 


s* 


a 


[ (ts*+r-62 )®+ 4 r (61+62 )1 

— — : ) 


+p( (1+26 i“ 62 ) (ts®+r-62 )+{2+2r ) (61+62 ) 


+ -Si (l+26i-6a )*1 . 
s 


(4.73) 


#* 


Due to ( 4 . 66 ), (4.67) and (4.37), we have 0 < eg 
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Wo nov7 use (4,37), i.Q. Q<t«lt and axpand (4.73) to terms 
of order t. Thon (4.73) may be written 


4r(6x+68 )+(r-6a )® 


9.i(P)9s(P) “ 


1 


2t£ (1+261-62 )(r-«3)+(2+2r){6j+6a).«j(t)] 

(pS+ II. - I iii.i I. yi ii. p 


4r( 61+63 )+{r-6a)a 


t* (1+261-62 )s 
4r(6t+6a)+(r-68)® 


and the zeros are easily found to be 


Pi* Pa» Pa 


r+6a i/ r6* 

Pi = — 5“ + 2i +o(t)} 

!5 8 

• 

-t[ (l+26i-6a ) (r-63 ) + (2+2r) (61+63 )-tO (t)l 

4r ( 61+63 )+(r-68 )® 


(4.74) 

(4.75) 


(-fii-fia ) {r*-4r6i+2E36j-26i63(l-r)+6|+6|+2ir6a+t»6i+0(t)) 

4- 2 ti ■■ ■ ■ ■■'■■■■■ 

4r(6|+6a)+(r-6a)a 

and pt is the ccn jugate ot Pt" Here we note that the 
expressions under the square root signs are positive due 
to (4.37). 


B. Determination of Branches for gi(p) and g2(p) 

Ail the zeros of gi(p)g2(p) are the zeros of 92 (p) 
cause a due to (4.63), (4.70) to (4.72), (4.74) and (4.75) » 


For computational purposes it is desirable to have p, , p^ 
in a series expansion of t. See Appendix 1. 
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we find that 

R® 9i(Pi) ” ) + o(t) (4.76) 

Iro gl(P2) > I* BtfTcx > 0, (4*77) 

Since gi(p) never vanishes on the entire p-plane with 
negative real -axis being deleted due to (4.70) to (4*72}« 
9i(p) not have branch points. Moreover# we find ftom 

equations (4.73)# (4.74) and (4.75) 

9a(p) - g'^fp ) ' (P*’Pi)^^^(P-Pi)^'*^^(P*Pa)^^^(p*Pa)^'^*» 

(4.7S) 

and this shows that pi, p^, p 2 are branch points of 
order <)ne for 9a(p). We define the domain of the 9 i(p) to 
be the entire p-plane with negative-real axis being deleted# 
and choose the branch gi(p) by the requirement that 

9i(p) = ^^(is 2-62 ) (;^+t)T(^^r+ty^^^ */ 2 r’ 

for p cs ^ 0 (4.79) 

We define a domain in which 9a(p) single valued such 
that the domain is the entire p-plane cut along the follow- 
ing curves ; 

(a) the negative real axis 

(b) the line joining Px and p^# that is# the locus of 
points p such that Rc(p) =* Ro(px) an4 

im(pt) < lm(p) :ilin(Pi) 

(c) the curve joining and p 2 # that iS/ the locus 
of points p along 
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( 1 ) the line such that Re(p) » Re(pa) and 
Ira(pa) < lm(p) < Re(pa)< 

( 2 ) the portion of circular arc whose radius is 

equal to « *713 Re(pa) and Whose argument lies 
between and ~ , 

(3) the line such that Re(p) » Re(pa) and Re(pa) ^ 
lm(p) < im(pS). 

We choose the single-valued branch of 9a (p) by the requite** 
ment that * 


9 a (p ) is 2 -6a )Yi+t ) + ( jgr +t V+2 1 63, - */lTr eTlw i+ 4 a 

for p = £ > 0 » (4 *80) 

Due tO^(4«79) and (4.80), we find 

i/a ‘ 1 

gi(p) * P9(l + *T77" P - 0 . GO (4.81 } 


92 (P) “ Ps(l “ + ^(|)) as p 

p '^s ^ 


00 


(4.88) 


We find that -€3 < Re(P 2 ) l«(Pa) ® with the aid of (4*37), 
(4.38) , (4.67 ) and (4 •75) as in the Figure 


Cat Formulation of w(^,p) in Convolution Form 

I 

It is expedient to define w^^*((,p), w^^^(C,p), 
and w^* such that 

w(C*p) = 


(4.83) 



• 92 



Figure 1. The p-*plane* 
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= s|r / w^*’^(?,p)exp(pT) dp with 1 »> 1«2 

(4.84) 


80 that hy the convolution theorem 

w(C#'t) “ / <^u* {4*85) 

0 . ’ . 

From (4«50)o (4*56)« (4.60)« (4*64) and (4.65) set 

W^‘^(C.P) = j 4 — {4.88) 

p 


wt*)(C,p) = I 


a, 


(p+b) 


,(65+62 -ps®+ 5 '{ (ps®-62 ) (p+t) 


(r-s®)p(p-a) ^ 

(p+h) ((ps®-6a )di+6id- } 


+ (pr'4*t)*{'2t6i ) } + 
Ai 


P9X92 


2(r-s^ )p(p**2i) 


Jexp (gi-Kta)C) 

( » .b 1 


-( 


Ai 


{6^ +62 -ps2 +|.( {ps» -62 ) (p+t ) 


(r-s^ )p(p-*a) ^ 

(p+t») {(ps®-6a)di+6^d2) 


+(ps:+t)+2t6i)'}+ 


Pgi9a 


Ax 

2(r-s2 )p(p-a) 


)exp{- (91- 92)0 




exp ( "*p^ /»C) (4.87) 

(r-s3)p(p-a) 

with a and b in (4.88) and (4.37) defined from (4.59) such 
that 

1 / 

Aa(p) T (s;“S*)pp '2(p-a}(p+b). 


(4.88) 
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Here a and '•b are zeros of 


3 ' 

ps*-6i-’62-t«pr-(pft)(ps3«a2)-2t62^«fpt+ptr « 0 
and they are 

a s — {r+ts*-s2-aa-t-tr-[ (s2-r-ts3+a2+t+tr )* 
2 (r- 3 ®) 

- 4 (r-s* ) ( -6® -62 -t+t6a -2t6i ) ) ^' 2 ) 


( 4 . 89 ) 


b a {r-f t 82 - 32 -da -t~tr-<-[ (s^ -t-te^-t-ea-t-t-t'tr )2 

'- 4 (r-s 3 )(- 6 i- 62 -t+t 6 a- 2 tdi)]^^ 3 ) (4 • 90 ) 

and we find that a > 0 , b > 0 by the relation ( 4 * 37 ). 

The inversion of w^^^(C,p) of ( 4 . 86 ) is found with 
the aid of the table [ 37 } 


N 


exp(-b(T-z^ )) dz. ( 4 . 91 ) 

To obtain a real integral representation of w^*^(^[#t) 
for T > Cf consider the integral 


2 




/ w^^(CiP)exp(pr) dp ( 4 . 92 ) 

evaluated along the closed contour shown in Figure 1 . 

D. inversion of wt®^(^,p) by Contour Integration 

\ 

V?e express w^^^(?,p) in a simpler form from ( 4 * 87 ); 


-(a)/- % -/ AaCp) 

w' '(s#p) =( — 


(r-s*)p(p-a) 9 x 92 2 (r-s* )p(p-a) 


^ )exp(-|;r!(9i+g3)C) 


•( 


A3(p) 


Ai 


(r-s 2 )p(p-a) ^ 1^2 2 {r-s 2 )p(p-a) 


Ax 

(r-s 2 )p(p-a) 


exp(«*p^^®5 ) 


)exp (- |^(gx " 92 ) C) 

^ X 


( 4 . 93 ) 
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whore 

AaCp) |<3i32(r-sa)p3 + (s3(|r«i+dia2« |6xd2 ) 

-f* 3:(6xda~ 1^162- |dxr)+ts2(. IdxS^'J^x'i^xr-tlaa ) )p2 
•«* |<326x)+ |d 26 x 62 * 2^x^a+t(s* (gdx-da 

+ |dira3-dir-dir6i+6ida+r6xda-dx62 


-^1^2 ) (^ 3 - 1 - 26 x )+ 2 (da-dx ) ( 6 x+ 6 a ) ) 

+ (da^-dx ) { '*6a*M+26x ) • (4.94) 

.With the aid o£ (4.81) and (4.82) we find that 
1/2 

exp(- fe— (gi+^a)? +pt) = 0(exp(p(T-O) as p 00 , 

” “ 1 . 




r^/zp^A 


(4.95) 


exp(- (gi-ga)C +Pt) = 0(exp(pT- £— ) ) 


as p — > 00 


(4.96) 


and consequently as p — > oo 

W^* ^(C.p)»8<p(pt) = 0(|exp(p(T-C))+ 0(|exp(pT- 


r 




)) 


+ 0(^Xp(pT-p^^*C))- 


(4.97) 


From (4.97) \m conclude that the contribution to the integral 
(4 .92 )N from the large circular portion of the contour goes^ 
to zero as the radius tends to infinity when t > C* From 
(4.73) and (4.93) we see that the contributions from the 
small circles around the branch points of ga(p) go to zero 
as the radii tend to zero. Let us consider the centribut’.on 



to t^a intC‘jraX from the small circle around the origin aa 
the rai^iiaa tend to zero* We may express of 

(4*0^) as 


w^^^(C.p) 


A3(p) 




(r-s®)p(p-a) ^1^8 


1 t % 

-r-IexpC- 




. - expf- §5 — (gi-ga)Ol 

" ^ 1 




2(r-sa )p{p-a) 


[exp(- §5^(9x+ga)C} 

** * X 


.V 


+ exp(- £|~(gi-gj )0-2exp(*p‘/®0]. (4.SS) 


We note that 


®^P(\ 2'f ^ ^ ^ ^ ^ •exp ( - 2£^ " (9i*9a )C) 


•. £ 


? — 92C+0{p) as p ■-> 0, 

ilL 

2fx 


(4.99) 


^ /s j / 

exp(- fs — (gi+ga)C)+«xp(- fy— (gi-ga )C)-2exp(-p '*C) 


s 2p^'»C(l- )+ 0 (p) as p ->• 0, (4.100) 

combining (4.98 ), (4.99) and (4*100), the contribution to 
the integral (4*92) along the small circle around the origin 
goes to zero as the radius tends to zero* Also, from (4.93) 
we see that the contribution from the small circle around 
P ** -’^2 9 oes to zero as the radius tends to zero because 
the integrand is bounded around p = -£2 * These considera* 
tions, with the aid of (4.89), (4*93) and Cauchy *s integral 
theorem , lead to 
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j-i"! (I^i / e>:p(pT)dp) 

s ^(C,t)+ ^ litn(/ + /+ /+/+/ +/ +/ 

ON NM ML KJ HGi GiG, FgPi 

+ / ■*•/+/+/ +/ +/ +/ 

PlE DC BA N'O* M'N* L'M' J'K' 


+ / + / + 

'IS 'II 

GiH’ G 2 G 3 , 

^ ^ 3 ( 3 ) Ai(a) 

(r-s 2 )agi(a)gj (a) 2a(r-s3) 

. Aj (a ) ^ Ai ( a ) 

(r-s* )agi(a)g 3 (a) 2a(r-s*) 

+ — *• ^ ^ exp(- a^^*5+ax) 
a(r-s 2 ) 


/ + / + / + / ) 

FiFj E'Fi CD* A'B' 

3 ^/ 2 - 

)exp(- 2 47 a J ( g 1 ( ^ ) ( ^ ) ) ) 

3 I /2 


(4.101) 


in which the integrand of the integrals in the brackets is 
^ (C*P)®^p(P'T' ) “lim*' refer to the limit process such 

that the large radius tends to infinity and the small radii 

1 / 

tend to zero. The values of p fj(p), 9i(p)/ 92 (p)# 
which are needed to evaluate the integrals in the parenthesis 
are to be determined consistently with the construction of 
the Riemann sheet described in (4-70), (4.71), (4.72), 

(4.79) and (4.80). 


D~l. Evaluation of gj^(p) along the Contour 

For this purpose it is expedient to introduce new 
functions 


Zi(p) = 9i(p) 


i = 1.2 


(4.102) 




(4.103) 


where wa define a cut for g(z. ) to be the negative real 
axis on the “plane such that 

G(i) = for £ > 0 

and (4.103a) 

G(“i) =* isTi for i > 0 . 

Then we shall choose proper signs along the contour in Pig. 1 
in the expression 

g^(p) « ± G(Z^) 


so that 9^(p) i = 1,2 are consistent with the Riemann 
sheet described in (4.70) to (4.72), (4.79) and (4.30), 

Due to (4.79) ’and (4.80) we may utilize' the Schwarz reflection 
principle, and in this case we have 


gi(P'^)' == (4.104) 

and since the contour in Fig. 1 is symmetric with respect 
to the real line, we determine the value of g^(p) along 
the contour which lies in the upper half-plane only and 
utilize (4.104) for the lower half-plane. Considering the 
mapping of the contour in the half-plane of Pig . 1 into the 
'1^-plane, and then into the G-plane under the restriction 
of (4.37), we find that 

9i{p) = G(Zi(p)) along HG 1 G 2 ,, F 2 F 1 B, DC, BA, JK, U\ (4.10$) 
gi(p) ~ “ g(Zi(p)) along NO (4.106) 

92 (p) ~ ^(^ 2 (p)) along HGj^G 2 » BA, iTK, £jM (4.107) 

92 (P) ~ “ g(Z 2 (p)) along F 2 p 3 ^E, DC, NO. (4.108 ) 
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We are now ready to evaluate the integrals in the parenthesis 
of (4.101). 

# 

D-2. integration along 0N-N*0“ for T i C and t < C 


Along ON, p « with £% < Ji < co ^ and by (4.X06) 
we find that 


^ ij ( 4 . 109 ) 

where 

X = - (r+ts^-62 )f+t+2t62-t62 

y = 2 r^/*s//*(i- 62 )V 2 (^_ej^)V 2 . 

The contribution to the integral of (4.92) along the 
contour ON-N’O' is then 


.® , -Ai(-£) 


1 / 


(r-s2 )B(i+a) 


sin (■/To + ( 


A3(-£) 


Ai(-^) 


)sin( 


./T 


2(r-s2 )i(i+a) s/TT7 

for T i. C and t < C* 


(r-s* )i(i+a) jgi 1 
(im gi)0) exp(-iT) di (4.110) 


D-3. 


Integration along the Arc NM-M*N* for t i C and 
X < C 


We next consider the contribution to the integral of 
(4.92) along the arcs NM~M*N‘ . Near p = we have 

1 / 

that p ^(gi-ga) 

For T < C* see 


Section D-7 . 
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are analytic due to the 


As (p) 


Ai(p) 


(r-s* )p(p-a)gigg 2(r-s3 )p(p-a) 
fact that along HL, p ® and 




-s2jj3-f{r+tsa-6g )jg-(t4-2t 6 i-t 63 )+ 2 a^/r , ( 4 . 111 ) 


93 = 9l 


and along OH/ 

9 a ® 9 x where g^^ is given by (4.109). 

Relations (4.52), (4.53), (4.54), (4.66) and (4.67) lead to 


/ 




or 


9i - 9a 4r /^2(p + eg) '2 


9i + 92 2r^^2p^/2(p 


(4.112) 


2fi 9i - 92 

and this last relation shows that the contribution to the 
integral (4.92) from the summand of the integrand, i.e., 

Va 

» 

1' 


As{p) 


(r-s® )p(p-a)gig 2 2 (r-s® )p(p-a) 

exp(~p^/»C+PT)/ 


■ )exp(“ |j~(9x+g3 )C+PT) 


(r-s® )p(p-a) 

vanishes as the radius of the arc NM*M*N* approaches to zero 
But the contribution to the integral (4.92) from the summand 


of the integrand, i.e., 

^ A 3 (p) 


Ai(p) 


p^'2 


(r-s®)p(p-a) 9 ig 2 2 (r-s® )p(p-a) 

is not readily established because fi(p) has a branch point 


at p = We consider the mappings 

u = (p € 1 ) 


(4.113) 
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which maps the contour in Fig« 2 onto the semi^circle C* 
in Tig* 3 with the direction clockwise in both cases. We 
del.ne 




^ / l(u) exp(- + (u* - 6 i)t) du (4.114) 


whe •« 


l(u) = 


2u 

2uAg (u»-€x) 


(r-s3 )(u®-6i )(u*-ei-a)gx(u»-ei)ga(u*-€i) 

UAi(u*-Ci) 


(4.115) 


(r-s*)(ua-€x)(ua-6x-a) 

v(u) = (u* - ei)^/a(gi(u®-ex)-gj(u»-6i)). (4.116) 

Kow we consider the branch of (p such that 

♦ - 

(P+€i)^^ 2 - - •/i+€i for p a: i > -Si, (4ell7) 


with the cut from *€i to -oo along the negative real axis* 
We consider the mapping 

u = (p + ei)^'*. (4.118) 

If we were to go over the circular arc twice in Pig. 2« 
then we see that the mapping (4*11S) will map the circular 
ate into the semi^circle C" shown in Fig* 3. As we did in 
(4*114)# we define 


" 57T ^ ^l(u)exp(- + (u»-ex)T) du. (4.119) 

c 

If we add (4*114) and (4.119)# then by the residue theorem 

« Residue (l(u)exp(- + {u*-€i)t)) (4.120) 

where u » 0 is the only singularity enclosed by the contours 
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C*-c”. Consider Now as given in (4*119) where the 

branch defined by (4.116) is used. Duo to (4.118)i (4ill9)f 
we see that the integration of (4.119) vanishes as 

the radius of the arc NM«M*N* approaches zero. 

p'plane u-plano 



Fig. 2 



D-4. Integration along KJ-or*K' for t i. C 


Along KJ# p ■ -f with 0 < *< 62 # (4.105) we 

find that 

g, . J X -K ^ (4.12X) 

where 

X » s*i* -(r+tB*- 62 )i+t+ 2 t 6 x**t 62 
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For X L the contribution to the integral of (4.92) along 
the contour KJ'*J*K* Is found to.be 


7 



n)0(~ 


Ai(-i) _ Aa(-i) 

r-8»)i(i+a) (r-s»)i(i+a)lgxj 


3 


) 


A f ^ 

- ain(»/TO ) exp (- tr ) d £ (4.122) 

(r-sa)l(i+a) 

where Ai(p), A3(p) and "a" are given by (4.58), (4.94) and 
(4.89). 


If 

D-5. Integration along ML*L*M* for t i 5 and t < C 


Along MLr p * -f with and with the aid of 

(4.105) and (4. ,107} we find that the contribution to the 
integral of (4.92) along the countour ML^Ii'H* is 



Ai(-je) 
(s*-r )^(i+a) 


ain(VI C) 


(4.123) 


for T C T < C * 


Integration along DC*-BA-A*B*-C’D* and HGi-<jiG 2"“ 
PaPl-FiE-E*PJ-riFi-^iGl-GiHl fox T i t 

Por T utilization of (4.104)^ (4.105) and 

(4.108) lead that the contribution to the integral (4.92) 
along the contour DC-BA-A*B*-C*D* and HGi-G3LG2-PaFi-PiE- 
ETi-PiFi-GiGl-ciHi vanishes. (4.124) 

For r < Ci s®® this section D-7. 
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D«>7* Modifications on integrations for t < C 


When T < the contribution to the integral of (4.63} 
from the summand of the integrand^ 


r ' ' ' 

(r-o» Jp(p-a)giga 


Al(p) 

2(r-s3 )p(p-a) 


)exp(- 



{«i+9a)C+PT) 

(4.12S) 


along the Bromwich contour vanishes because we let the 
contour for this summand be closed by the right arc of the 
circle which is shown as a dotted curve in Fig* 1* Then 
the oantribution from this circular arc to the integral of 
(4el25) approaches to zero as the radius gets large. Hence 
for r < some modifications should be made on (4.122) 
and (4«124)« but the results in (4*110)« (4*120)« (4.123) 
ere not affected. 


D-8. Integration along DC-BA-A’B* and 

PaPi-PiE-E'Pi-PlPa-GjGi-GiHi for T < c 


Along DC« p - Re px with 0 ^ f ^ Xm Px 

an4 we let 

l/o 

p ' » * ax + 

(p + €3 ) ® bx + ibj (4.126) 

(p + €x)^^^ * Cx + ic^. 

' /fith the aid of (4.52)« (4.53)* (4.102) and (4*103)« we have 
Zx • Xx + iYx 


where 


(4.127) 
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Xi •" s*(Re*Pi-jg* )+r+te *-62 )ro Pi+t+2t6i«t6a 

+ 2r*^/*8(ci(aibx-a2ba )-C2(aiba+a2bi)) 

Yi = 2is*Re Pi+^(r+ts*-62)+2r^/2e^c^^aj^b3-aabi) 

and 

Z 2 » X 2 ^ iYa (4.128) 

where 

Xa « 8 * (Re*Px-^* )+(r+ts*- 62 )Ro Pi+t+2t6i-t62 
“ 2r ^*8(ci (a^bx-aaba )“C 2 (a^ba aab^)) 

Ya r 2is*Re Pi+i(r+ts»-6a )-2r^^®s{Cx(axba-aabx ) 

+ Ca(aibx-aaba)). 

Frem (4 •105) and (4«103)t we have that 
9x = g(Zi ) 

and the contribution to the integral of (4.92) along DC»BA 
-A*B*-C*D’ when *r < is 


1 I« Pi 
7 Jo 


exp{- giC+pT) x/ 

Re — > ■ ■■ ’ ■» ■■■>(Ax (p)sinh (^~ gaC) 


(s®-r)p(p-a) 

- cosh g 2 C))]d£. 

9x'=ia 


(4.128) 


■Along Fa?!, p • Re pj ij with im P2 ^ t ^ 2^'>Re pj , and 


by (4.105) end (4.108) 


9i =G(Zi) 
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9a a= “<5(2a ) 

where Zi, are given by (4*127) and 4*128)i and the con« 
tribution to the integral of (4*92) along 
for T < C is 




j ^3Re Pa exp(- ^ giC+Px) _ 

/ Re[ — -» . T ■ ^ ■■ ■■■(Ax(p)ainh(g'g^ gaC) 

im Pa (s*-r)p(p-a) 




2A3(p) 

9i92 


o^/a 

cosh(~ gaC))] di * 


(4.130) 


Along GiH» p « paexp(i0) with 0^9^ 37 t/ 4, and 

by (4.105) and (4.107) 

9i • G(Zi ) 

9i - G ( Za ) 

where 2^ , Za given by (4.127) and 4.128)/ and the 
cohtribution to the integral of (4.92.) along HGi-P|.E*E*Fl 
-GiH* for T < C is 


,i/a 


2"/*(Re Ps) o exp(- § 5 ^ giC+pT) i/, 

/ . (sine[im( E"Ai(p)sinh( 2 -^--g 8 C) 


sir/. 


(r-s» )p(p-a) 


p^^2 


■ ^g®gg' 92^)1 }1 


+ COS 9 (Re { 


p^/s 

exp(- giC+pt) 

(r-s2)p(p-a) 


lA 


(Ai(p)ainh(gY^ gaC) 


2 A 3 ( p ) /s 

^ -5T3T- 


gi9a 


(4.131) 
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D-9. Integration along KJ*J'K* for r < C 


The contribution to the integral of (4.92) along KJ- 
J*K' for T C ia easily deduced from (4.122} as 


i ® A. ( 

r / «5ip(*‘iT)sin(s/TO ~ — — 

^ €a (r-s^ )f (jB+a) 


d^. 


(4.132) 


D-10. Obtained by inversion of 

We consider two cases, i.e., t ^ ^ and t For 

T i C# have that from (4.101) 

fa}/« - 1 rr A3(a) Ax(a)^ , ^ 

^u.T) - HrTrT(9x(a) 

, •> , As(a) Ai(a)^ , , , v 

+ ga(a))C«T) “(g 7( - a 7 g T- ( a)' + 

- 92 (a) )C+aT)+Ai(a)exp(-a’'/®5+aT)l - Ii(CiT) (4.133) 

where 1 i(C,t) = ^ lim( /+/+/+/ 

ON NM ML KJ 

+ J* + f ■ +/ **'/ ) 

K*0* M*N* 1,‘M* J*K* 

whose integrals are evaluated in (4.110)i, (4.120), (4.122) 
and (4.123). For t < C# we have that from (4.101) and the 
vanishing of the integration of (4.125) 


,(0 


■*■ )exp(- ^|^j 5 j(gi(a)”g 2 (a) )C+aT))-l 2 (C.t) (4.134) 

Where 
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Ja(t.T) 


1 


2wi 


lim( / 
ON 


/ 'J* / + / 

m ML KJ HGi 


^ I + / 

GiGs PaPi 


+ / *^/*^/'^/ '*■/ 
FlE DC BA N’O* M*N 


+ / ^ / + / 

L*M* J»K* GiH* 


+ / 9+ / + / ) 

GaGl FiFa S*f1 C»D' A*B* 

'hose integrals are evaluated in (4.110), (4.120), (4.123), 
4.129), (4.132}i (4.120) and (4.131). 


E. Inversion of w(^,p) in Real Integral Form 

Combining (4.133), (4.134), (4.S5) and (4.86) leads 

:o the complete description of the displacement field of 

;he solid component w(?,t) 

^ 0 

O T -/u \ 

/ exp(-b(u-z*'* ) ) da '(C»t-u) du. (4.135) 

*^^0 0 ■ 

proro (4.133), (4.134) and (4.135) we see that w(5it) 
satisfies the boundary and regularity conditions specified 
by (4.33) and (4.34). 

Since the material constants have to be determined by ex** 
periment for the mixture and such an experiment has not 
yet been devised, we ^all not attempt any further invest! •• 
gation about the behavior of the displacement field of the 
solid component at thi^ point, even though we have the 
exact solution given by (4.135) which may be evaluated 
numerically by computers. 
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4,5. Ea^Xy Time Solutions 

The numerical evaluation of (4,135) does not seem to 

'be an easy task. One way to avoid this difficulty is to 

represent v;(5,p) in a power series with respect to — for 

P 

sufficiently large p, and then invert the resulting expres 
Sion term by tern* This procedure leads to an early 

tine solution for . 

As p CD , we have that 


2r^/asp^/2(p+£i)^'^a 

• = ^ !!S. * * «1 * ik + o(Ml 


P 


P‘ 


(4.136) 


Where 


Ml 


Ms 


3 ^ (t+tr)-r( 6 | 6 a) 


2 rs' 


(* 1 + 62 )i 2 rs® (t+tr) ( 61 + 62 )+s^ (t+tr )* ) 
8r®s^ 


M3 


I , C 6 !+£a)® (6!+63)(t+tr)' 




r®s^ 



'M4 


- l>5(t+tr)^ 20(t+tr ( 6 t + 60 ) 

+ — — — - 


2{t«tr)2(6?+6a)8 


2 k 4 


r^s 


5 (t+tr ) (61+62 ) 

' "i* ««<* 


rs 


s 


e 
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„ 1 , 7(wr)8 . 5(t^tr)M^i +<5a) 2 (t+tr )3 ( 6 ^+^ )» 

“» " 256^ x» ■ r*s3 “ • x30* 

2 (t+tr )® {fij+ig )® 5(t+tr)(6*+ij)^ ‘^( 61 + 62 )®% 

+ . ■" '' ■_ - ., , «, ~v • 

jc^sO rs® 

With the aid o£ (4.68), (4.69), (4.79), (4.80) and (4.136), 
we have that as p — > o3 


where 


Ni 
N 2 « 

N 3 =* 

N 4 = 

Na » 


ps[l 

Ni 

■** 

P ® 

52 . ^ 

p 

N3 1 

pp/2 



Ns 

^ Nio 


+ + 

9 * 

p 4 p 

+ — — 
/8 p8 

ne 1 1 

Wl 

5 i 

N3 

PS 

A 

P 

" ” ' i 7 a 

pp '® 


+ 2s. - 

N9 

1 

p-p 

AmUSL 

/fl -.6 

*r 

3 

r 

'■/a 




N 5 


P® P®P 




^7 

p4p"' 


+ + o(~)j (4.137a) 

D« 


N4 


^6 Ng 

27 T *:7 


Ny 


P* ' P®P 


— + 0^~ 
p®p '2 p® 


p-p 


1/3 


3r>*'a 


(ts* - 62 ) 



2s^ 

t(l+26i-62) 


k: 


2s® 


- NjNa - 




- N1N4 - NaNa 



m 




He '■•NiNg - N 3 N 4 - 




N, 


s 


- NiHa - NgHg - NaN^ 


He « -NiN 7 - NaNa - NaNg - 


2 


Ns 




- NiKa - N 2 N 7 - Natle “ NaNg 


nS 


Nxo ® “N 1 N 9 - NgNs - N 3 N 7 - N 4 N 5 - 


r^/ 2 j 4 . 


Nn * 


- NiNio - NaNft - N3N8 “ N4N7 - NftN^. 


As p 00 , w© have that | ^ p"y iS| due to (4<54)> (4«66}« 

^ (i... 1 


1 6i . 3 el 3‘5 €j 

3 5 T 5 ^ . 

P® p- 


^ITpT " 77i;^":“ 2 -7 * 8 :7 " 5T^ „3 


3’5'7 £4 3. 5. 7. 9 e| 


4 '. 2 '* p 4 


5 ‘. 2 ® 


~ + 0(— )). (4.138) 

a t^e 


P“ P 

Equations (4.137) and (4.138) lead to, as p -<> oo , 

exp(- ff^*(gi+gj )0 = expl-pC-(Na- |^)tHl“ ^(Ne* j^iHa+fcv) 


+ —^[“N8 + 
P® 

N4C 1 

3Nae| 

3-5 3 ^ 


2 

8 

+ ' €1 + 

3123 

2-(N4- -5— 

+ ^[-Ne+ 

Ne€i 

3N4e| 

3*5 3 

3-5*7 4 

2 

3 

•*• — ■ Nae,' 
3123 

4 '.24 


8 


^iNa 32 ., 3 2 3*5 3 ^ 

+ (N4 -f §^i>(Ng + gNjCj- —■"■ C l )C 

* ■* 3 *. 2 ^ 


(Nj-EiNa-*- i- € 1 )* 

TT 


C *1 + 0 (^ )]. 


(4.139) 
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.‘/a 


exp(- (gi-ga) C)*expl-p — — - 


(Hs-exNi/2)C 


(Hs-CiWi/2)*ca [(fsiej- H5 )C+(Hs- 5«i«i)*C»1 
* 3p 


(**#“ 5**i*i ) (§**i*i“ §N*ej. + Ng)5‘ 




Al(p) 


l«t US consider v ' / ' p^ given in (4.56)> (4 .57} and (4. 58} 
As p —> GO we have 


sJtIt ‘ “ ^(Ra-Rs-RiRa) 


pap‘/2 


. i,{2BaR3-Ra>aRS-RxRa). ^(R^3K?R,.R^28,H.8,.HaR,» ) 

f 

t- ^(•38983 •HR|R3<Ra<*>RaR|-3RiR^R, 48183) 4 o(~)] (4.141} 


P 

Where 


Ri 


(dj (t+6j. )/di}*6^'-t 

8*-r 


Rji == s*-r+62+t(l+r-s® ) 

“ *5^*62**t(l*62+263L ) . 

From (4*73) we immediately have 

gi{p)5a(P) « P*s«(l + £i. + S|. + S|. + £i) 

P P* P* P^ 


(4.142} 


2(ts*-r-62 ) 


8 


a 


where 
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Ca = 


2t(-l+25i-6a-2r) (tsa+r-fij )*+4r ( 6^+62 ) 


s* 


s' 


^3 26j,** 53 ) {tfi^+r "*62 )+(2+2r ) (6^ 


C4 = ~(1 ■*“ 261-62 

With the aid of (4 •142)^ ve find the asymptotic expansion 

57(pfe‘(p) 

giTp)9a'(p) ^ ^ 

• * |(C^2C,C3) 


i|c!ca ^O(^)). 

^ 4!2* P 


‘(4.143) 


With the aid of (4.143), we find the asymptotic expansion 
of b(p)‘ and D(p)# given in equation (4.50), is 

di Si S: 

P' 

where 


B(p) 


1/1 ®i ®® 1.. 


(4.144) 


St = 


6id2-di53 (ts2+r-6a-2s* ) Cj, 


djsa 


2 s* 


•C2 32 Cl 2(6id2-di6s) (tsa+r-6a-2sa ) 

s* = - -^ + ^1- r-< — *^i+*^* ^2— — 

(t(l+ 25 i- 6 j )+2 (61+62 ) + (Ri-R2)(ts*+r-a2"2s») 


2s® 
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Sa 


{-Rj+Ra ) (t (1+261-62 )+2 ( 61+62 ) ) 

2 e 8 

(RI+K 3 +R 1 R 2 ){ts 2 - 6 a+r- 2 s* ) Ci o 

+ • - ^{-Ra+Ra+RiRa 


2 s^ 


( t ( 1+2 61 —62 )+2 (61+62 ) ) (Ri "R a ) ( ts® —62 +r— 2s® ) 


8 


a< 


/“Ca 3 a . 2 ( 6 id 2 -^i 62 ) (tB®+r- 63 - 2 s® ) 

+ (- 4 - + ^ Ci)( Ri+Ra- 




a 


+ (“ 4 ^ + fCxCj *• I? Cl)» 

d(p) \r- (Do •*■ “ + “§■ + o(~)) 

2p3p^'a Pi*® 


(4. US) 


where 


2(6id2-di6a) (ts®+r- 62 - 2 s® ) Ci 

Do • 2 R 1 - 2 R 2 + + 2 "* 


diS® 


C 2 3 3 Cl 2 ( 6 id 2 -di 62 ) 


Di *= 2Ra-2R3-2RiR3+ 5 g- Ci + 


2 


dis' 


- R 1 +R 3 


(ts®+r- 62 - 2 s® ) 1 2 

2 )+ — 5“( t ( 1 +2 6 1 -6a ) +2 ( 61+63 ) 


8 


+ (Ri-Ra )(ts®-6a+r-2s» )) 



US 


Dj => 2(2RaR3-R|+RiRj-RiR3)+(Ri-Ra} 


(t(l+26i-6a)+2(6?+6a)) 


(-Ra+Ra+riRa)(ts®+r- 6 a- 2 s 3 ) Cj , a 
( t (l+26i-6a )+2 ( fix+^8 ) ) ) (ts®+r-6a -2s* ) 


,3a Ca _2(aida-di6a ) (tsa+r-6a-2s* ) 

" (sCi- S-)(” ; Bi+Ra 

d|^s® 

/ *"3 3 _ _ 5 „3 . 

- (“ 5 “ + 4 CiCa - ^ Cl ) . 

Contbining (4.50), (4.56), (4.60), (4.139), (4^140), 
(4.141), (4.144) and (4.145) and inverting term by term, 
we have that for early time 

w(?,t) “ di(-4T)®'^3i®erfc +(Ra-Ri ) (4 t )*^ 3i®er5e 
+(R 3 +RiRa-Ra)( 4 T)’/® i^erfc ^ + 

3 5 

■»exp((-Sa'*‘ |^)OH(T”O( f( 02 T '~ 

.5(N,-^.|e“)).I^f^(Sa-SiC(l»4-^ I ■ 

./ N 4 C 1 3 . 2 5 3 (N4-€iN 2/2+ |- e|)^ 

+;(-Ns+ ^ N 2 €i + ~ €1 + 


0 + ‘ 

NiC 


§^( 4 t)^^® i®erfc — — 4 - " N 3 )( 4 t)* i®erfc 

^ aA ^ ^ 2 A 


(Di+ ^(Ng- ^ NiC 

' " * " ' * ' (4t) i^erfc — — + •• 

. 2/V 


(4.146) 
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Hero orfc(x) is the compllmentaj^ error fUnistibn defined 

toy 


CD 


er£c(x) ^ ) 


and the repeated integrals of the conplemiantary error fune* 
tien are defined by 

i^erfc(x) *» erfc(x)» 


i%rfc(x) “ / i"“** erfc(t) dt# n » 

X 

See [38] for example* 

A similar procedure may be used to find an early time 
solution for v(^,r). We begin by finding the asymptotic 
expansions of .the factors in equations (4.51)# (4*57) as 
p 00 • ^The following factor is written as a series 
expansion in terms of (^) 

(p(^-^p)(ps*-^8)+t6i)- |p(gi+g2)») K Ka ' i . , , 

2 Ko+ — + -y 0(~) (4.147) 

(p(ts*+r6t )+t( 6^-62 ) ) P P P 


00 


n-i 


where 


Ka 




t( 6 i -62 )-s® (2N4+Na ) 
ts*+r 6 i 

-2s»(N8+NaN4) t(6i“6a )(t(4i“6a )“8®(8M4+JI* ) ) 


Ka = 


ts®-fr 6 i 

t 3 ( 5 i- 62 )a 


(tfi**fr 6 i )‘ 


(ts»+t4i) 


, 2 taa( 6 i- 6 a) 

a (t(«i“ 6 a)- 8 *( 2 N 4 +Na))+ — ~T~(** 6 +Na» 4 ) 


(ts*+r 6 i ) 


- a^( 2 Ne + 2N2N6 + nJ) 
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The following factor is written as a series expansion in 
terms of (i) 

Jr 


(p{(t-t-p)(ps»-63)+tai)- |p(gi-ga)*) 
(pts*+r )+t( 6 i -“62 ) ) 

where 


l>oP*+I<xP+I>a+0(|) 

(4.148) 


Lo 


Ll 


ts^+rdi 

ts*-62-s2Nx ts^{6x-6a) 


ts^+r6x 


(ts2+r6x 


t2s*(6x-62 )® t(6i-6a)(ts»-62-saNi) t(6i-63 )-2 s3NiN3 

Lj + — — 


(ts*+r6x )® 


(ts^+rfii ) 


(tsa+rfii) 


the following factor is written as a series expansion in 
terms of (^) 

« ‘ Ai fo) (ps^“6a]d, 6xda 

?(4i+6a-ps®+t6i+(p8*-6a }(p+t) + ly^' 


p(ts*+rfix )+t(6x-6a ) 

— 171/ "a r(Jo+ ^ ^ 

pp 'a(tsa+r6i,) P P* P® 


(4.149) 


here 


3 » 6* (Rx-Ra )+tS®-s2- ^ 


X ®* s* (Ra ■•Ra^'RxRa “Ra )+ts^ •62 "•s® ) + (^x'^^a +t6x*^^a ) 

d, t(6i-6a) 

- (s»{Bi-Ra)+tB»-8»- g- 6j)- 


(ts^+r6x ) 



ns 


t® (61—62)^ do ■ 

“ 7 — ^(sa(Ri-Rs)+tsa- 8 a- ^ i^) 

(ts®+r6i) 

) ‘ 2 ■ ■ 

(a^ (Ra-Rg-RiRa )->(Ri-Ra ) {ts*-6a-'«* ) 

(ts2+r6i) 

+ ( “^.^2 ) )‘^a^ (^RfiRs^Ra ^RiRa *Rx^a J 

+ (ts2-a3-s3) (R|.R3-.RiRJ)•^(R^-Ra)(6J+6att(6l-63). . 

With. the aid of (4.51), (4.57), (4.61), (4.139), (4.140), 
(4.141), (4.145), (4.147), (4.148) and (4.149)> and with 
the elementary inversion process, we have that for early 
time 


v(C,t) « di[ 


7 J 0 


(t^2+roi ) 


(4^)^/a iierfc- 


2^ 


Ja 


« ^ (4T)^^^i^erfc — ^ - 

(ts^+rdi) (ts^+r6j,) 

+exp {(-H2+ 1 ) C }h ( t-c ) t ~' r ~ ( 5 7 2'} — + " " r ' (772 y~ 

T2(t-C)^^ 2 Zg .l/„ : , NiC 

k" "'* +1 + ^(4 t) i*^erfc 

2/T 


’(4t)*^^ i*er£c^ 
2^ 


rW^ 


2^ 3/ ^ NiC 

I i®erfc 






+ •-*] , 


(4.150) 


where 

.To == I<o 

Ti ® Ki + Ro(Si**C(N4“ I* eiNg I* € 1 )) 


-»! I 
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“ Kj+Ki t g^tNa* §■ ®i) )J'^ptSa"SxC(tl4“ j€jiNa+ ^i) 


+(-Hj+ 58461 - |na£i)+ ^ Ci+ 


(N 4 - ItJaei + § ct)* 




Zo = D«Lo . 

Zi - - C^o^( 83~ ^ £lNi), 

Z 3 * .I-ofDi + 50 oC*( 83“ sNiSi)*) + DoLi. 


We note that these early time solutions o£ w(C>t) and 
v(C»t) are in effect without any restrictions beyond the 
conditions of (2.40) and (2.41). And the early time 
solutions of all other field variables follow trivially 
from the equations (4.27)« (4.146) and (4 .150). 



SUmARY A£n> CONCZAJSZOHS 


Xn this thssio# we have reviewed the major oontrihu* ’ 

tipns to the development of a theory of meohanically and 

thermally interacting continuous media. Beginning with the 

work of Darcy and Terzaghi we have traced the work of Biot# 

* 

xruesdell and Toupin and the recent work of Green. Naghdi# 
Steel. Atkins and Chadwick. As is usual, the theeretical 
development has preceded the number of applications and in 
thi'i thesis we have al&empted to utilize a linearized 
version of the mixture theory to derive results which are 
readily applicable to practical boundary value problems. 

Our first result is in the form of an integral rela» 
tion commonly known as a reciprocal theorem. Xt relates 
the solution of one problem to that of another problem 
each of which is due to different boundary and initial 
data, we have indicated how this theorem reduces to a 
theorem applicable to a single constituent and we have 
shown how one might use such a theorem, indeed, we intend 
to explore its uses in future research much along thcf lines 
vaed in classical elasticity. 

Our second result cansists of a solution of a 

fundamental initial boundary value problem using the 
linearized mixture theory. It is the first actual boundary 
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value problem to be solvsd using mixture theory* Due to 
its connlexity the results are given in integral form only* 
Further development must await experimental evidence con- 
cerning the size of the material properties* Such experi- 
ments, incidentally f are a second possible line of future 
research and it is our intention to attempt to devise 
simple analytical models which will lead to estimates of 

the material constants* We will be guided by those methods 

« 

used in single constituent theories* 

The integral representation of the solution of the 
boundary value problem given in Chapter 4 is exact .to 
terms of order t* but, due to the complexity of the inte- 
grands in the integrals, not much can be Inferred about 
the displacement field* For this reason we have given the 
starting solution, i«a*, the ea^ly time approximation* 

This solution may prove more useful as far as actual e<xnpu- 
tation is concerned* 
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APPENDIX 


The location of the zeros of (4.73) may be given in 
power series of t as follows, 
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Pi = {^) (^ - ^a) + 0(t3)) 
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and th^ir conjugates pt and pj » where w© used the fol 
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